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^  In  a  laser  triggered  spark  gap  (or  spark  column)  breakdown  is 
initiated  by  a  laser  which  creates  a  narrow preionized  column  between  the 
electrodes.  This  small  plasma  column,  of  diameter  *•  50  Atm,  expands  to  an 
arc  with  a  diameter  of  0. 5-1.0  irm  on  a  time  scale  of  100  ns,  and  is  capable 
of  carrying  a  current  in  excess  of  10  kA.  During  the  expansion  phase,  the 
resistance  of  the  plasma  column  decreases  from  many  kilo-ohms  to  as  little 
as  a  few  tenths  of  an  ohm.  For  many  applications,  the  duration  of  the 
current  pulse  through  the  spark  gap  is  only  100-200  ns  and  therefore  is 
cotranensurate  with  the  time  required  for  the  plasma  column  to  expand. 
Therefore  during  the  time  of  interest,  the  spark  column  changes  from  being 
a  dominant  resistive  loss  to  being  a  lossless  circuit  element  in  comparison 
to  the  load.  In  order  to  optimize  the  performance,  and  minimize  the 
losses,  of  laser  triggered  spark  gaps  operating  under  these  conditions,  it 
is  desireable  to  have  a  detailed  understanding  of  the  physxcs  of  the 
expansion  phase  of  the  spark,  ftiis  report  susmtarizes  the  first  phase  of 
development  of  a  model  to  achieve  that  goal.  We  discuss  a  model  which, 
from  first  principles,  describes  the  growth  of  the  arc  in  a  laser  triggered 

spark  gap  in  a  gas  mixture  of  a  noble  gas  and  a  diatomic  molecule  ,  (*-  - 

(eg.  i  Ar/H2).  The  model  makes  no  assumptions  as  to  thermal  equilibrium  of 
the  plasma  column.  Heavy  particle  hydrodynamics,  electron  kinetics, 
radiation  transport,  and  an  external  circuit  are  included  in  the  analysis. 
Skin  depth  affects  for  the  penetration  of  the  applied  electric  field  into 
the  plasma  are  also  included.  Results  from  the  model  are  discussed,  tfe 
find  that  the  expansion  of  the  spark  column  is  primarily  by  convection  from 
the  hot  ionized  core.  An  ionized  channel  of  low  mass  density  is  confined 
within  a  high  mass  density  compression  (or  shock)  wave.  The  convective 
expansion  of  the  core  is  augmented  by  thermal-  and  photo-ionization  of  the 
neutral  gas  at  the  interior  of  the  high  mass  density  shell.  Based  on  the 
results  of  the  model,  and  a  thermodynamic  analysis,  a  simple  thermodynamic 
model  for  a  laser  triggered  spark  gap  is  presented. 
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SECTION  I 


INTRODUCTION 


Spark  gaps  are  traditionally  triggered  by  overvoltage  or  by  internal 
auxiliary  electrodes.  The  use  of  lasers  as  trigger  sources  for  spark  gaps 
has  significant  advantages  over  these  methods* 1  *  Laser  triggering  of 

spark  gaps  has  been  demonstrated  as  a  method  whereby  jitter  can  be 
minimized.  For  applications  where  many  spark  gaps  must  be  triggered 
simultaneously,  lac  ?r  triggering  provides  a  method  whereby  this  can  be 
accomplished. 

In  laser  triggered  spark  gaps,  a  laser  is  used  to  initiate  the 
electron  avalanche  and  breakdown  of  the  gas.  A  laser  with  sufficient 
intensity  to  ionize  the  gas  is  focused  between  the  electrodes  of  the 
device.  (See  Pig.  1-1)  Although  the  potential  across  the  spark  gap  may  be 
insufficient  to  cause  self  breakdown,  once  a  plasma  column  is  initiated  by 
the  laser,  the  voltage  is  sufficient  to  sustain  and  expand  the  column.  The 
duration,  location  and  intensity  of  the  laser  can  be  carefully  controlled, 
and  therefore  so  can  the  time  of  breakdown  ox  the  spark  gap. 

Although  many  systematic  studies  of  laser  triggered  spark  gaps  have 
been  performed*1'1-1'  *  and  some  modeling  has  teen  done* 1 ' 4* 1 ' 1  ,  a  first 

principles  model  of  a  laser  triggered  spark  column  has  not  teen  presented. 

A  first  principles  model  would  treat  all  the  basic  phenomenon  without 
making  any  apriori  assumptions  as  to  the  thermodynamic  ntate  of  the  gas, 
and  would  include  hydrodynamics,  electron  and  heavy  particle  kinetics, 
radiation  transport  and  an  external  circuit.  Development  of  such  a  model 
is  described  in  this  report.  We  have  attempted  to  fully  describe  all  the 
pertinent  physics  in  the  plasma  column  without  prejudice  to  a  particular 
preconceived  notion  of  the  dominant  processes.  Although  this  method  is  not 
the  most  computationally  expedient  technique,  it  is  the  most  revealing. 

Reference  is  made  in  this  report  to  experimental  measurements .  These 
measurements  were  made  under  contract  to  the  Naval  Surface  Weapons  Center 
(Contract  N60921-83-C-A057)  concurrent  to  the  development  of  this  model. 

We  acknowledge  the  support  and  advica  of  NSWC  in  the  experimental  program. 

The  model  is  described  in  detail  in  Sec.  II.  Simulated  spark  column 
characteristics  are  discussed  in  Sec.  III.  (We  use  tho  phrases  "spark  gap" 
and  "spark  column"  synonymously. )  The  rates  of  arc  formation  and  expansion 
are  discussed  in  Sec.  IV.  The  simulation  of  interferograma  from  results  of 
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the  model,  to  compare  with  experimental  data,  is  discussed  in  Sec.  V. 

Based  on  the  results  of  the  experimental  and  theoretical  studies,  a  simple 
thermodynamic  model  for  a  laser  triggered  spark  gap  is  described  in 
Sec.  VI. 


Figure  1-1.  Geometry  for  Laser  Triggered  Spark  Column 
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SECTION  II 


DESCRIPTION  Or  THE  MODEL 


The  model  •xplicitly  treats  a  spark  column  with  an  initial  gas 
mixture  of  argon  and  molecular  hydrogen.  The  energy  levels  Included  in  the 
model,  discussed  below,  were  chosen  accordingly  for  these  species.  Sparks 
in  other  gas  mixtures  though,  can  be  approximated  by  scaling  the  energy 
levels  and  atomic  masses  used  in  the  model  from  those  of  argon  and  hydrogen 
to  those  of  the  desired  species.  To  approximate  a  diatomic  molecule  other 
than  hydrogen,  the  mass  of  the  desired  molecule  (and  atomic  const ituient ) 
is  used  in  place  of  that  of  hydrogen  where— ever  mass  appears  in  the  model. 
The  energy  levels  of  molecular  hydrogen,  and  their  excitation  rates,  remain 
unchanged.  The  energy  levels  of  the  atomic  constituient  are  scaled  in  the 
same  fashion  as  that  for  substituting  a  different  rare  gas  for  argon. 

To  approximate  a  rare  gas  other  than  argon,  the  same  mass 
substitution  as  described  above  is  used.  The  appropriate  electron  momentum 
transfer  cross  section,  thermal  conductivity,  and  ionization  potentials  for 
the  new  gas  are  also  used.  The  energy  levels  appearing  in  the  code, 
selected  for  argon,  are  scaled  in  the  following  manner.  The  threshold 
energy  of  the  first  excited  state  is  revised  to  that  of  the  new  rare  gas. 
The  ionization  potential  of  the  first  excited  state  is  then  used  as  the 
scaling  parameter.  The  separation  in  energy  between  adjacent  levels  is 
increased  or  decreased  in  the  ratio  of  the  ionization  potential  of  the 
first  excited  state  for  the  desired  rare  gas  as  compared  to  that  of  argon. 
The  revised  values  of  the  energy  levels  are  used  in  the  calculation  of 
excitation  rates  as  required . 

II .A.  Heavy  Particle  Collisions 

The  importance  of  heavy  partical  collisions  which  result  in  the 
Changr.  in  the  level  of  excitation  of  one  of  the  collision  partners  is 
particularly  important  with  respect  to  expansion  of  the  arc.  The  arc  can 
expand  by  only  one  of  three  mechanisms.  A  partially  ionized  fluid  element 
can  convect  from  an  inner  radial  point  to  an  outer  radial  point,  electrons 
can  diffuse  from  an  inner  radial  point  to  an  outer  radial  point,  or  local 
sources  can  provide  the  seed  electrons  required  for  avalanching  and 
subsequent  rapid  ionization  of  the  gas .  Two  major  local  sources  exist  for 
these  seed  electrons;  photoionization  and  thermal  ionization. 


As  the  gas  temperature  is  raised  in  the  core  of  the  plasma  column, 
the  increase  in  thermodynamic  pressure  provides  the  driving  force  for  an 
expansion  wave  to  propagate  radially  outwards.  At  the  outer  limit  of  the 
expansion  wave  is  a  shock  or  compression  wave,  ttie  compression  wave  is  a 
thin  region  of  gas  with  a  nass  density  that  is  large  as  compared  to  the 
quiescent  gas  beyond  the  shock.  To  a  good  approximation,  this  gas 
undergoes  an  adiabatic  compression  and  therefore  has  an  associated  increase 
in  gas  temperature.  For  a  Ap/p^  (ie.t  change  in  mass  density  divided  by 
the  initial  gas  density)  in  the  compression  wave  larger  than  3,  the  gas 
temperature  is  increased  to  a  value  in  excess  of  1000  K.  At  these 
temperatures  thermal  ionization  of  the  gas  (as  opposed  to  direct  electron 
impact  ionization)  is  no  longer  a  process  that  can  be  ignored.  Although 
the  rate  of  thermal  ionization  is  many  orders  of  magnitude  less  than  that 
by  electron  impact  ionization,  its  importance  is  manifested  by  the  fact 
that  in  the  absence  of  free  electrons  ( and  ignoring  for  the  moment 
photoionization),  thermal  ionization  is  the  only  local  source  of  electrons 
in  this  region.  Once  local  free  electrons  are  generated  by  thermal 
ionization,  and  in  the  presence  of  a  sufficiently  large  applied  electric 
field,  the  more  rapid  rate  of  electron  impact  ionization  dominates  that  of 
the  thermal  ionization  and  initiates  a  local  electron  avalanche .  The  non— 
convective  expansion  of  the  plasma  column  progresses  radially  outward  as 
the  electron  density  of  the  newly  avalanched  region  approaches  that  of  the 
core  of  the  plasma.  The  electron  avalanche  also  serves  to  rapidly  heat  the 
heavy  particles,  thereby  raising  the  local  thermodyisaic  pressure  and 
perpetuating  the  convecting  expansion  wave. 

Concurrent  to  thermal  ionization  in  the  adiabatically  compressed 
shock,  thermal  dissociation  and  vibrational  excitation  also  occur. 

Locally,  thermal  dissociation  and  vibrational  excitation  are  not  important 
in  directly  providing  an  impetus  for  the  electron  avalanche  to  occur.  The 
importance  of  thermal  dissociation  is  manifested  indirectly  by  chances  in 
electron  transport  properties  that  may  occur  as  a  result  of  the  relative 
change  in  atomic  and  molecular  densities.  For  example,  electron  impact 
excitation  of  vibrational  modes  of  the  molecular  species  is  one  of  the  more 
efficient  processes  by  which  electrons  lose  energy  and  hence  by  which  the 
electron  temperature  is  constrained  to  a  relatively  low  value.  If  in  the 
adiabatically  heated  region  a  large  fraction  of  the  molecular  species  is 
thermally  dissociated,  the  rate  at  which  electrons  lose  energy  as  a  result 
of  vibrational  excitation  will  be  reduced.  This  reduction  enables  a  higher 
average  electron  temperature,  a  larger  rat#  of  electron  impact  ionization, 
and  therefore  a  more  rapid  avalanche  of  the  gas.  Analogously,  the 
importance  to  the  avalanche  process  of  thermal  vibrational  excitation  of 
the  molecular  species  is  not  direct;  rather  it  is  an  indirect  effect 
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manifested  by  changes  in  the  rate  of  electron  impact  ionisation.  The 
ionization  potential  of  a  v irrationally  excited  species  is  smaller  than 
that  of  the  ground  state  molecule  by  a  value  equal  to  that  of  the 
vibrational  excitation  energy.  For  a  given  electron  temperature,  the  rate 
of  ionization  of  the  vibrationally  excited  molecule  is  therefore  larger 
than  that  of  the  ground  state  molecule.  A  large  fraction  of  molecules 
vibrationally  excited  by  thermal  collisions  can  accelerate  the  avalanche 
process  by  providing  a  reservoir  of  molecules  with  a  lower  ionization 
potential  than  that  of  the  ground  state  molecules.  Although  not  direct 
contributors  to  the  avalanche  process,  thermal  dissociation  and  excitation 
of  the  molecular  species  can  indirectly  be  responsible  for  a  rapid  rate  of 
electron  avalanche  by  the  influence  the  products  of  these  processes  have  on 
the  rate  of  electron  impact  ionization. 

The  rate  of  thermal  dissociation  of  molecular  species  in  the 
adiabatically  heated  gas  was  obtained  from  the  literature  and  is  given  in 
Table  11-2.  the  rate  of  thermal  ionization  was  approximated  as  being  gas 
kinetic  with  an  Arrenhius  factor  exp(-B^/kT  ),  where  Ej^  is  the  ionization 

potential  of  the  heavy  particle  and  T  is  tfie  local  gals  temperature. 
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ZI.B.  Photoionization  and  Radiation  Transport  ! 

Photoionization,  in  analogy  to  thermal  ionizatio[i,  is  not,  on  an 
absolute  scale,  an  important  source  of  ionization.  Thu  importance  of 
photoionizatlcn  is  that  it  provides  a  local  source  of  ionization  in  the 
neutral  gas  exterior  to  the  plasma  column.  As  discussed  above,  local 
sources  of  ionization,  independent  of  local  electron  collision  processes, 
are  necessary  to  provide  the  seed  electrons  required  for  an  electron 
avalanche.  An  electron  avalanche  initiated  by  local  siources  of  ionization 
in  once  neutral  gas  is  the  only  non-convective  mechanism  by  which  the  arc 
may  increase  its  radius,  and  is  the  only  mechanism  by  which  new  mass  can  be 
entrained  into  the  expanding  core. 

The  radiation  which  produces  photoionized  seed  electrons  in  the 
region  exterior  to  the  arc  originates  in  the  plasma  within  the  arc.  The 
flux  of  radiation  at  location  r  of  frequency  v,  <t>(  r ,  v ) ,  resulting  from 
volumetric  sources  of  radiation  at  location  r’,  S(r',v),  is  given  by 
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where  a(z'',v)  is  the  photon  absorption  coefficient  for  frequt.icy  v  at 
location  r*'.  The  integral  in  Eg.  IZ-1  (called  the  photon  transport 
integral)  has,  in  general,  no  analytic  solution  if  either  the  source 
function  or  absorption  coefficient  are  complex  functions  of  radius.  For 
our  conditions,  S(r,v)  and  a(r,v)  are  not  functions  of  a  single  variable 
( eg . :  temperature )  but  rather  represent  the  sum  of  many  processes  which 
depend  on  the  electron  density  and  temperature,  and  the  density  of  a  large 
number  of  excited  states.  The  integral  in  Eq.  IZ-1  Must  therefore  be 
performed  numerically  for  an  "exact"  solution.  For  conditions  where  either 
S(r,v)  or  a(r,v)  change  rapidly  for  the  time  scale  of  interest,  repeated 
numerical  evaluation  of  Eq.  II-l  is  costly  and  approximate  solutions  to 
Eq.  II— 1  must  be  used.  Approximate  solutions  can  be  obtained  provided  that 
S(r,v)  and  a(r,v)  are  constant  over  some  finite  region,  or  the  integral  in 
the  exponent  in  Eq.  IZ-1  is  neglibly  small.  Both  of  these  criteria  hold  to 
some  extent  for  our  conditions  and  are  exploited  to  obtain  approximate 
solutions . 

To  a  good  approximation,  photons  generated  within  the  plasma  column 
can  be  classified  as  being  optically  thick  or  optically  thin.  Optically 
thick  photons  are  those  photons  that  have' a  high  probability  of  being 
absorbed  over  the  path  lengths  of  interest,  for  our  conditions,  these  are 
photons  of  sufficiently  high  energy  that  they  are  capable  of  ionizing 
ground  state  species.  Optically  thin  photons  are  those  photons  having  a 
low  probability  of  being  absorbed  over  the  path  lengths  of  interest.  For 
our  conditions,  these  are  photons  of  low  enough  energy  that  they  are 
incapable  of  ionizing  ground  state  species,  but  are  capable  of  ionizing 
excited  states.  These  two  photon  groups  naturally  arise  from  the  form  of 
the  photoionization  cross  section  and  the  energy  separation  of  excited 
states . 

For  hydrogenic  levels,  the  phot  ionization  cross  section  scales  as 
(€/hv)  ,  * 1  ‘ 1  ^  where  €  is  the  ionization  potential  of  the  level  and  hv  is 
the  photon  energy.  Therefore  photons  with  energy  far  in  excess  of 
threshold  have  a  small  photoionization  cross  section  for  that  level. 

Recall  that  for  hydrogen  and  most  noble  gases,  the  ratio  of  the  ionization 
potential  for  the  ground  state  as  compared  to  the  first  excited  state  is 
about  3  to  5.  Therefore  photons  which  efficiently  ionize  the  ground  state 
do  not  have  a  large  cross  section  for  ionizing  even  the  first  excited 
state.  Photons  that  are  near  threshold  for  ionizing  excited  states,  and 
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therefore  have  a  large  cross  section  for  doing  so,  are  not  energetic  enough 
to  icnize  ground  states.  By  grouping  photons  into  energy  groups  which 
satisfy  the  optically  thin  and  optically  thick  criteria  described  above, 
two  different  approximate  solutions  of  Eq.  II-l  can  be  applied. 


For  optically  thin  photons,  the  exponential  function  in  Eq.  XX-1  can 
be  ignored  and  the  integral  reduces  to 
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where  ©■  (r )  is  the  photon  flux  at  location  r  for  optically  thin  photons 
grouped  into  an  energy  range  denoted  by  k.  We  used  cylindrical  coordinates, 
assumed  azimuthal  symmetry,  approximated  that  the  axial  dimension  is  large 
compared  with  the  radial  dimension  and  that  the  photon  sources  were 
constant  over  in  a  cylinder  of  thickness  Ar  and  radius  r.  With  these 
conditions. 
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where  the  summation  in  Eq.  I 1-3  is  over  cylindrical  shells.  The  function 
f(t)  is 
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For  optically  thick  photons,  we  approximate  that  a(r,v)  is  constant  over 
tl  e  path  between  the  point  of  emission  and  the  point  of  observation.  For 
the  same  geometry  the  optically  thin  case,  we  define  the  parameters 
a^-r^/R,  b^Ar^/R,  and  c-z/R.  The  integral  in  Eq.  II-l  reduces  to 
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where  the  average  absorption  coefficient  of  the  path  length  of  interest  is 
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The  upper  limit  o£  integration  for  the  axial  direction  is  chosen  such  that 
(i/R)/8  *  1.  The  integral  in  Eq.  II-5  was  performed  numerically  and 
parameterized  as  function  of  a,  and  $.  The  results  were  compiled  into  a 
lookup  table  (function  g  in  Eq.  II-5)  for  use  during  execution  of  the 
program. 

Por  our  plasmas  of  interest,  the  absorption  coefficient  is  the  sum  of 
two  dominant  terms.  The  first,  and  of  lesser  importance  in  the  partially 
ionized  phase,  is  inverse  bremsstrahlung  or  free-free  absorption.  The 
aborption  coeffiecient  for  photons  of  frequency  v  in  a  plasma  with 
temperature  T  is^II#1* 


<*„(  v )  -  E 

i 


2  6 


.  , _ ,  n  N.Z. e 

4  f  2n  1  e  i  i 

3  l3kTj“  3/2.. 


3/2  3 
hem  v 
e 


*ff 


.1 


[11-7] 


where  the  factor  gff  *  1, 


H ..  is  the  density  of  ion  i,  and  is  the 
ionization  state  of  ion  i.  The  second  and  dominant  absorption  mechanism 
during  the  partially  ionized  phase  is  bound-free  electron  tranritions,  or 
photoionization.  Por  hydrogenic  states,  the  photoionization  absorption 
coefficient  for  frequency  v  is  approxmately 
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where  JT.j  is  the  density  of  species  i  excited  to  level  j  and  is  the 
ionization  potential  for  that  state.  In  practive,  photons  were  grouped 
into  energy  bins  with  widths  of  a  few  electron  volts.  When  doing  so  the 
absorption  coefficient  for  a  particular  energy  bin  was  the  average  value  of 
of  Eq.s  II— 7  and  II-8  for  the  range  of  energy  of  interest.  When  the  plasma 
becomes  fully  ionized,  and  at  longer  wavelengths,  inverse  bremsstrahlung 
becomes  the  more  important  absorption  mechanism. 


I.C.  Radiation  Source  Terms 

Although  line  radiation  is  included  a.?  an  energy  loss  mechanism 

(ie.:  spontaneous  emission),  radiation  transport  in  calculated  only  for  the 

continuum  component  of  the  spectrum.  The  radiation  source  function  S(r,v) 
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(photons/ cm  -s)  for  this  continuum  component  results  from  free-bound 
radiative  recombination  and  brems3trahlung ,  The  former  source  is  the  more 
important  source  at  lower  electron  temperatures  and  densities.  The  latter 
is  the  more  important  source  at  higher  electron  temperatures  and 
densities.  Assuming  hydrogenic  levels,  the  source  function  for  free  bound 
radiation  resulting  from  recombination  of  electrons  with  temperature  Tg 
recombining  into  a  level  with  ionization  potentional  €q  and  emitting 
radiation  with  energy  between  and 
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where  the  summation  is  over  all  atomic  and  molecular  states  in  the  plasma 
which  satisfy  the  energy  criteria  c.  >  e  y  »  m _/(2€  . ),  and  0  « 
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m  /(2kT  ).  N  is  the  density  of  atoms  in  the  next  higher  ionization  state 
them  the  level  to  which  recombination  is  occuring  and  A(Te)  is  a 
normalization  constant  to  insure  that  the  sum  of  photon  sources  equals  the 
rate  of  radiative  recombination. 
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The  source  function  for  free-free  (bremstrahllung)  radiation 
resulting  from  electrons  with  temperature  T0  colliding  with  ions  of  charge 
state  z  and  density  N2  and  emitting  radiation  with  energy  between  and 
«...  i.*11-2’  1  % 
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in  Eqs.  Il-ix  to  11-14,  r(Z,Te,  e^,  c^+1 )  i&  the  rate  constant  for  the 
process  by  Maxwellian  electrons  with  temperature  Tg  and  CT^^(Z,  £ ,  6^,  e^+1 


is  the  analogous  cross  section  for  a  single  electron  having  energy  ( . 
Eq.  11-12  is  the  averaging  of  the  cross  section  for  energy  {  over  the 
Maxwellian  electron  distribution  function. 
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II.D.  Hydrodynamic  and  Thermodynamic  Equations 

1 

*;  In  this  section,  we  will  describe  the  hydrodynamic  and  thermodynamic 

equations  used  to  model  the  plasma  column.  These  equations  treat  the 
plasma  as  a  single  fluid  but  with  two  temperatures .  The  nearly  fully 
ionized  core  in  the  center  of  the  arc  is  close  to  being  in  thermal 
equilibrium  and  a  single  temperature  might  suffice.  Exterior  to  the  core, 
the  electrons  and  heavy  particles  are  not  in  thermal  equilibrium  and 
therefore  two  temperatures  are  required,  the  electron  temperature  and  the 
heavy  particle,  or  gas,  temperature.  Conceptually,  each  temperature  is 
associated  with  a  given  mass  fraction.  In  similar  situations,  one  would 
require  two  fluids  to  accurately  describe  the  pertinent  hydrodynamics, 
electrons  and  heavy  particles.  For  our  conditions,  though,  the  two 
particles  are  sufficiently  well  coupled  that  a  single  fluid  convective 
velocity  suffices.  This  fact  results  from  the  high  ion  density,  low  rate 
of  diffusion,  and  high  heavy  particle  collision  frequency. 

For  our  conditions,  the  ion  density  is  sufficiently  high  that 
electron  transport  is  accurately  described  by  the  ambipolar  approximation. 
Using  this  approximation  (see  below  for  details),  the  diffusive  motion  of 
electrons  is  constrained  to  that  of  the  ions,  the  rate  of  which  is 
increased  a3  a  result.  The  region  of  charge  separation,  though,  is  small 
compared  to  our  spatial  resolution.  Therefore  the  density  of  electrons  and 
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ions  are  locally  equal.  The  low  rate  of  diffusion  for  all  particles,  as 
compared  to  the  convective  motion,  insures  that  the  ions  convert  at  the 
same  velocity  as  other  heavy  particles,  the  high  collision  frequency 
between  heavy  particles  insures  that  all  heavy  particles  can  be  described 
with  a  single  temperature. 


In  the  anibipolar  diffusion  approximation, 
can  be  written  as'11,1^ 


the  radial  electric  field 
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where  Pg  is  the  thermodynamic  electron  pressure  and  ng  is  the  electron 
density.  The  force  of  this  electric  field  acts  only  on  the  ions.  However, 
due  to  the  high  collision  frequency  between  heavy  particles,  the  momentum 
acquired  by  tho  ions  is  distributed  among  both  ions  and  neutrals  in  a 
nearly  uniform  fashion.  The  contribution  of  the  racial  electric  field  to 
*  ’"omentum  equation  is  therefore  weighted  by  the  mass  fraction  of  ions 
ied  by  their  charge  state. 

With  a  single  fluid,  albiet  with  multiple  temperatures ,  there  is  no 
relative  motion  between  species.  The  convective  velocity,  though,  is 
sufficiently  large  that  .'he  linear  dimension  of  the  region  of  interest  may 
change  by  more  them  an  order  of  magnitude  during  the  discharge  pulse  While 
the  entrained  mass  changes  by  a  significantly  smaller  value,  tinder  these 
conditions,  a  Lagrangian  formulation  of  the  transport  equations  is 
appropriate.  Using  this  formulation,  the  fluid  mass  within  each 
computational  cell  remains  constant  While  the  cell  boundaries  move  with  the 
convective  velocity.  In  this  manner,  the  density  of  mesh  points  tends  to 
increase  in  regions  of  compression  Where  gradients  are  liXely  to  be  large. 

The  fluid  continuity  and  momentum  conservation  equations  used  in  the 
model  are 
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In  Eq.  11-17,  the  sum  is  over  heavy  species  and  is  the  charge  state  of 
species  i.  The  last  term  in  Eq.  11-17  is  due  to  the  effect  of  self 
generated  electromagnetic  fields  and  is  the  source  of  the  compressive  force 
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in  "pinch"  plasmas .  For  our  collisionally  dominated  conditions,  this 
effect  is  not  important. 

Separate  energy  equations  were  written  for  the  electrons  and  the 
heavy  particles.  Energy  is  coupled  into  the  plasma  through  the 
acceleration  of  electron  by  the  applied  electric  field.  Heavy  particles 
are  heated  by  collisions  with  electrons.  The  energy  equations  are  coupled 
to  the  momentum  equation  through  the  pressure  gradient. 
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The  electron  energy  conservation  equation  is 
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The  single  subscript  refers  to  a  particular  state  of  a  given  species  and 

the  double  subscript  r. .  for  a  rate  constant  or  energy  refers  to 

transitions  or  energy  differences  between  species  i  and  j.  In  Eq.  II-18,  Pg 

is  the  partial  pressure  of  the  electron  gas,  is  the  electron  thermal 

conductivity,  J  is  the  current  density,  O  is  the  electrical  conductivity, 

is  the  energy  separation  between  species  i  and  j,  and  e*  is  the 

ionization  potential  of  species  i.  The  rate  constant  for  excitation  of 

to  N,  is  r. .,  the  rate  constant  for  ionization  of  N.  is  r1,  the  rate 
1  i]  1  j+. 

constant  for  collisional  radiative  recombination  of  ion  N.  to  form  species 

/  )4>  C  1 

N'  is  r. .  and  the  analogous  rate  constant  for  radiative  recombination 

1  x  T 

is  r^.  is  the  electron  collision  frequency,  and  o^v)  is 
photoionization  cross  section  of  species  i  for  frequency  v. 


The  tenti3  in  Eq.  11-18  are  for  heating  of  the  electrons  by  adiabatic 
compression  of  the  electron  component  of  the  fluid,  thermal  conduction. 
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joule  heating,  thennalization  of  electrons  with  the  heavy  particles, 
inelastic  and  superelastic  excitation  collisions  with  heavy  p&rticles, 
ionization  collisons  with  heavy  particles,  electron  heating  due  to 
collisional  radiative  recombination,  energy  loss  due  to  radiative 
recombination,  energy  loss  due  to  bremsstrahlung  radiation,  and  energy 
inputs  to  the  distribution  resulting  from  jahotoionization .  In  the  term  for 
bremsstrahlung  loss  the  constant  fc  «■  10-1  eV  1  ^  2-crn  /s  which  is  small 
compared  to  joule  heating.  Tne  thermal  conduction  term  serves  to  smooth 
the  electron  temperature  distribution  as  a  function  of  radius,  however 
because  electron  transport  is  dominated  by  collisions  with  heavy  particles, 
this  term  is  not  a  dominant  contributor  to  the  total. 

Recombination  heating  is  the  transfer  of  thermal  and  potential  energy 
of  a  recombining  electron  to  a  second  electron  during  the  collision  which 
occurs  during  collisional  radiative  recombination.  Radiative  recombination 
differs  from  collisional  radiative  recombination  in  that  in  the  former 
process  the  kinetic  and  potential  energy  of  the  recombining  electron  is 
released  only  by  emitting  a  photon.  Recombination  heating  is  a 
particularly  important  term.  In  the  core  of  the  arc  where  quasi-steady 
state  conditions  are  obtained,  the  rate  of  recombination  is  equal  that  of 
ionization,  and  the  dominant  form  of  recombination  is  collisional 
radiative.  For  illustrative  purposes,  assume  that  the  second  colliding 
electron  carries  off  on  the  average  twice  the  thermal  energy.  Then  for  an 
electron  density  of  10  and  electron  temperature  of  4  ev,  the  rate  of 
recombination  heating  is  IS  GW/cm  ,  which  is  comparable  to  the  rate  joule 
heating,  approximately  90  GW/cm  .  Recombination  heating  is  not  a  true 
heat  source  in  the  sense  of  being  an  external  source  of  energy. 
Recombination  heating  returns  to  the  electron  distribution  the  kinetic 
energy  of  electrons  previously  "lost"  as  a  result  of  inelastic  ionizing 
collisions.  This  energy  is  stored  as  the  potential  energy  of  electrons  in 
the  continuum  relative  to  the  excitation  level  of  the  state  from  which  the 
ionization  oc cured . 


The  heavy  particle  energy  conservation  equation  is 
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where  cp  is  the  heavy  particle  specific  heat,  X  is  the  heavy  particle 
thermal  conductivity  and  is  the  electron  collision  frequency  with 
species  i.  The  terms  of  j  are  for  heating  from  adiabatic  compression,  heat 
transport  due  to  thermal  conduction,  and  joule  heating  by  collisions  with 
hot  electrons.  In  the  last  term  of  j,  for  changes  in  internal  energy  due 

n 

to  heavy  particle  collision  processes,  r^^  is  the  endothermic  rate  at 
Which  reactants  i  and  j  collide,  resulting  in  product  species  k  and  1.  The 
activation  energy  for  the  process  is  The  reverse  exothermic 

reaction  results  in  heating  of  the  gas.  Examples  of  these  heavy  particle 
collision  processes  are 


+  M<3>  -  H<*>  +  H(,°  +  M<J>, 


-  4.5  ev 
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H2(v  -  0)(i)  +  M<:^  -  H2(v  -  l)<jt)  4  M(1>,€i3kl  -  0.5  eV  [11-21] 

I 

where  M  is  any  second  or  third  body.  The  forward  heavy  particle  reactions 
are  especially  important  in  the  dissociation  of  the  molecular  species. 

This  results  from  the  rapid  heating  of  the  gas  by  electron  collisions  and 
the  relative  inefficiency  with  which  electrons  dissociate  molecules  as 
compared  to  heavy  particles  of  the  sar-.1  temperature . 

j 

II. E.  Kinetic  Equations 

The  hydrodynamic  and  thermodynamic  equations  described  in  the 
previous  section  yield  the  fluid  density,  heavy  particle  temperature ,  and 
electron  temperature  at  a  given  radial  point.  The  specific  properties 
required  for  solution  of  the  hydrodynamic  and  thermodynamic  equations 
(ie.:  energy  loss  rates,  electrical  conductivity,  thermal  conductivity)  are 
obtained  by  solution  of  kinetic  equations  which  provide  the  density  of 
individual  species  within  the  given  computational  cell.  From  these 
densities  the  necessary  thermodynamic  properties  and  reaction  rates  are 
obtained.  The  species  included  in  the  model,  as  functions  of  radial 
position  and  time,  are  listed  in  Table  II-l.  For  atomic  neutral  and  ionic 
species,  the  notation  (n)  in  Table  II-l  indicates  an  excited  state.  For 
hydrogen,  the  first  ten  states  of  atom  were  used.  Due  to  the  large  number 
of  excited  states  for  neutral  and  ionized  argon,  actual  states  were  grouped 
into  pusedo  states,  comprising  the  actual  states  lying  inside  an  energy 
interval  approximately  1500  cm  to  3000  cm  wide.  This  represents  a 
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somewhat  coarser  selection  of  states  thai  that  used  for  calculation  of 
partition  functions.  (See  Sec.  VI.)  The  statistical  weights  of  the  pseudo 
states  are  the  sum  of  the  statistical  weights  of  the  actual  levels  falling 
within  the  energy  interval.  The  energy  levels  used  for  Ar  I  to  Ar  VI  are 
listed  in  Table  II-3 .  Although  properties  are  included  in  the  model  for  up 
to  five  times  ionized  argon,  most  of  the  results  discussed  here  are  for 
calculations  where  the  highest  ionization  state  included  is  Ar  III.  The 
density  of  Ar  III  for  these  conditions  is  sufficiently  small  that  higher 
ionization  states  are  not  necessary. 

For  computation  of  the  electron  temperature  both  vibrational  and 
electronic  states  of  H2  were  included.  However  only  the  v  -  1  vibrational 
state  and  H2  (c  »  8.85  eV)  electronic  state  were  carried  as  kinetic 
variables.  For  our  collisionally  dominated  plasmas,  th>  lifetimes  of  these 
states  are  very  short  as  they  relax  to  the  ground  states  (or  first 
vibrational  state)  quickly.  Inclusion  of  the  electron  impact  process  for 
excitation  of  the  state  While  not  including  the  state  as  a  kinetic  variable 
is  equivalent  to  having  an  instantaneous  relaxation  of  the  state. 


The  processes  included  in  the  kinetics  portion  of  the  model  are 
summarized  in  Table  II-2.  Conservation  equations  for  each  species  were 
written.  The  entries  in  Table  II-2  are  for  processes  which  change  the 
density  of  a  particular  species  as  a  result  of  electron  or  heavy  particle 
collisions.  The  change  in  density  of  a  particular  species  due  to  non- 
collisional  processes  (ie.t  convection)  is  a  quantity  derived  from  the 
hydrodynamic  equation  for  heavy  particle  density.  In  the  Lagrangian 
formulation,  the  hydrodynamic  equations  provide  the  mass  density  in  each 
computational  cell,  while  maintaining  a  constant  total  mass  in  the  cell. 
Therefore  the  fractional  change  in  mass  density  in  the  cell,  resulting  from 
expansion  or  contraction  of  the  cell  boundaries,  manifiests  a  similar 
fractional  change  in  the  particle  density  of  each  species  within  the  cell. 
Therefore  a  correction  term  for  convective  effects  is  included  in  the 
conservation  equation  for  each  species.  A  typical  conservation  equation 
for  a  kinetic  species  in  the  model  is 


^  n.(r3i“j  '  rijHi]  -  neriHi  4nel<  *  rrne]H+ 

*  ".(‘KPil  4  8i„)  4  2  ra(K2  4H2<v))t,il  4  ra<(6li  4  6„t)] 
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In  Eq.  11-22,  0^  denotes  that  the  process  only  pertains  to  the  hydrogen 

excited  state  n  *  i  and  denotes  a  third  body  which  has  associated  with 

it  a  unique  rate  constant  also  having  subscript  k.  The  following  terms  in 

Eq.  11-22,  in  order  of  presentation,  have  rate  constants  that  are  function 

of  electron  temperature:  electron  excitation  of  from  level  (rate 

constant  r. .);  electron  collisional  relaxation  of  H.  to  H,  (r, electron 
31  “t  +.  i  j  -13 

impact  ionization  of  (r^  );  recombination  of  B  populating  by 

radiative  recombination  (rz)  and  by  collisional  radiative  recombination 

C  T  1  f 

(r.  );  dissociative  recombination  of  populating  H(l)  and  a  higher 
excited  state  (r*j);  electron  impact  dissociation  of  H2,  a2(v),  and  K2 
(rd>.  The  following  terms  in  Eq.  11-22,  in  order  of  presentation,  have 
rate  constants  that  are  functions  of  heavy  particle  temperature: 


collisional  quenching  by  heavy  particles  of  a  higher  lying  level  j  to 
(r^. )  and  of  H.  to  a  lower  level  j  (rh. );  thermal  dissociation  of  H  , 
«Jt>.  h"  and  H*  (rate  constant,  r£3.  r^<v».  r"  *.  r»*\,  and  2 


“2'  /  •  2  9  2  ' 

reassociation  of  H  atoms  to  make  H 


dk '  dk  'dk  'dk  - 

The  last  term  in  Eq.  11-22  is  the 


correction  to  the  density  of  resulting  from  convection. 


II. F.  Circuit  Equations 

From  a  system  point  of  view,  the  spark  column  is  merely  a  reactive 
element  in  a  discharge  circuit.  The  driving  force  for  the  kinetics  and 
hydrodynamics  for  the  spark  column  is  the  applied  electric  field  from  the 
discharge  circuit.  The  coupling  between  the  microscopic  properties  of  the 
spark  column  and  the  measure  able  macroscopic  parameters  is  described  by  the 
external  circuit  model. 

A  schematic  of  the  external  circuit  model  is  in  Pig.  II-l.  The 
capacitance  of  the  equivalent  external  circuit  is  CR,  its  inductance  l^, 
and  its  apparent  resistance  R^,.  The  spark  column  of  length  1  has 
resistance  Rg  and  inductance  Lg.  The  parallel  resistance  Rp  (•  looo  fl)  is 
that  associated  with  diagnostic  probes.  The  resistance  of  the  plasma 
column  is 
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t 


„  .<3/2  _  j(r) 

*  /  2vt^T)  at 


[ 12-23 J 


where  d  is  the  column  diameter,  j(r)  is  the  current  density,  o(r)  is  the 
electrical  conductivity  and  Z  is  the  total  current  through  the  column.  The 
column  resistance  is  equivalently  written  as 
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An  Vrl* 


2nrdr 


where  n@(r)  is  the  electron  density,  m^  is  the  electron  mass,  and  vc  is 
the  electron  collision  frequency.  The  spark  column  inductance  is 


.  «;‘/2  L_LL> 

10 


H  R 

o.  c 

*77  "  l(  d/2  ) 
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where  the  magnetic  field  at  location  r  is 

B(r)  *  rrJT  J(r* )  27/r‘dr’ 

0 


[ZI-26] 


In  Eq.  11-25,  Rc  is  the  radius  of  the  current  return  path.  The  total 
current  through  the  spark  column,  X,  is  the  integral  of  the  current  density 
j  over  the  cross  sectional  area  of  the  plasma.  The  local  current  density 
is  simply  written  as  ^  -  )E(  r )  where  M#(r)  is  the  electron  mobility  and 
E(r)  is  the  local  elec  r\,*  ield.  For  our  conditions  E(r)  is  truely  a 
local  variable  and  thexeiore  must  be  obtained  by  solution  of  Maxwell's 
equations.  (That  is,  E(r)  is  not  necessarily  V/£.)  The  electric,  or 
equivalently  the  magnetic  field,  applied  to  the  spark  column  must  "diffuse'* 
into  the  interior  of  the  plasma.  This  effect  is  known  as  skin  depth.  The 
local  current  density  is  therefore  obtained  from  solution  of 


mn .  u  n  +  £_i  ±  frJLfjLL£2] | 

dt  dr  4ffr  3r  iI3rta(r)JJ 


[11-27] 


c  1  3 
4JT  r  3r 


(r  ^(u(r)B(r))] 


With  the  quantities  defined  above,  the  equations  solved  for  the  external 
circuit  are 
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V  -  V 


at  c  '  at 

c 
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where  the  voltage  drop  across  the  plasma  column  is 


s.a* 

at 


V  1  V 

i  — £  +  *  «__e 
Rpj  *p[at  Ls 


[11-29] 


For  conditions  where  the  Lp  is  small,  -  I( RgRp  )/< Rs+Rp )  ■ 


II. G  Transport  Coefficients 


II.G.l.  Electrical  Conductivity  and  Electron  collision  Frequency 

The  electron  collision  frequency  enters  into  the  model  in  at  least 
three  places;  electrical  conductivity,  electron  thermal  conductivity,  and 
momentum  transfer  to  heavy  particles  (ie.i  gas  heating).  In  a  partially 
ionized  plasma,  the  electron  collision  frequency,  is  the  sum  of  the 
collision  frequency  with  neutral  particles,  v  ,  and  with  ions,  v  .  The 

Cn  Cl 

former  frequency,  is  given  simply  by  the  product  of  the  neutral 

particle  density  with  the  convolution  of  the  momentum  transfer  cross 
section  and  the  electron  distribution  function,  in  this  case  a  Maxwellian 
characterized,  by  the  temperature  T@ .  The  later  frequency,  vcJ,  is 
normally  given  by  the  Spitzer  value. 


'«  '  l  W 


477  c  m 
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where  the  sum  is  over  ions  in  *"he  plasma  having  Charge  and  density  N^. 
When  the  ion  collision  frequency  is  large  compared  to  the  neutral  collision 
frequency,  the  conductivity  approaches  the  Spitzer  limit.  <IX-3>  in  this 
limit,  the  conductivity  is  only  a  first  order!  function  of  electron 
temperature.  There  is  a  weak  dependence  on  electron  density  as  a  result  of 
the  appearance  of  n0  in  the  expression  for  A.  I 
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The  conductivity  of  interest  is  that  value  for  notion  perpendicular 

to  the  magnetic  filed,  a..  This  value,  as  opposed  to  o„  is  used  because 

■l  ■  /  tt  3  \ 

our  plasma  is  highly  collisional  and  satisfies  the  criteria1  ' 

13  3/2 

1.2  x  1C  BT  / 

Vvc - 2tnAn  — 1  ♦  1  t11-32’ 

e 

where  u>c  is  the  electron  cylcotron  frequency,  v c  is  the  electron  collision 
frequency,  B  is  in  Gauss,  T@  is  in  eV  and  n@  is  in  cm~  .  For  our  typical 
values, wc/vc  *  2  x  10-  * 


The  derivation  of  Eqs.  11-28  and  XX-29  assumes  that  there  are  many 

3 

electrons  within  a  Debye  sphere,  or  equivalently^  (4/3)jtXd  . »  1  *•  A  »  1. 

This  condition  requires  that  n  «  3  x  10  T  '  cm  ~  ,  where  the  electron 

mi  3)  e  e 

temperature  is  in  eV.  *  '  For  our  plasmas,  this  conditions  is  not 
necessarily  met.  A  first  order  correction  to  Eq.  XX-30  is  to  replace  inA, 
which  in  itself  is  an  approximation,  with  the  exact  expression. 


inA 


-J* 


0  (i*,3) 


—  dx>  — 
2  2 


in 


(l  +  A2]  --A 


1  +  A 
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Further  refinement  in  the  value  of  vcJ  requires  detailed  kinetic 
modeling.  For  non-Debye  plasmas,  the  conductivity  is  usually  lower  than 


the  value  one  would  calculate  in  the  Spitzer  limit, 
refinements  in  this  value  have  been  made. * 11 ,4* 


Semi-emperical 


I X . G . 2 .  Thermal  Conduct ivity 

Thermal  conductivity  enters  into  the  calculation  in  both  the  heav'’ 
particle  and  electron  energy  equations  and  is  a  mechanism  whereby  spatial 
variations  of  the  internal  energy  of  the  plasma  are  smoothed  in  the  absence 
of  convection.  The  conduction  of  heat  in  the  heavy  particle  fluid  is 
important  only  at  the  boundary  of  the  hot,  low  density  core  and  the 
relatively  cold,  high  density  shock  wave.  Conduction  of  heat  from  the  core 
into  the  cool  shock  contributes  to  thermal  ionization.  Recall  that  it  is 
only  by  thermal  ionization  or  photo-ionization  and  subsequent  electron 
avalanche  that  the  plasma  column  cam  expand  by  non-convective  means. 

The  electron  thermal  conductivity  is  a  function  of  the  electron 
density,  temperature,  and  collision  frequency  with  heavy  particles/11,1* 
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where  v  is  the  electron  collision  frequency  with  all  heavy  particles. 

Cn 


The  analogous  expression  for  heavy  particles  cannot  be  readily 
written.  Although  the  thermal  conductivity  of  a  high  temperture  gases  has 
been  measured,  we  require  the  value  for  arbitrary  mixtures  of  Ar  and  h2  in 
non-equilibrium  conditions.  We  therefore  must  resort  to  the  first 
principles  expression.  The  thermal  conductivity  for  a  mixture  of  heavy 
particles  is 


s  l  8  s 


[— r/a*  i  - 

IjtmJ  sJ  ] 


97  -5 
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2M  1  1/2Q  | 

sr  j  *sr 
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where  Ng  is  the  density  of  species  s,  Mg  is  the  mass  on  species  s,  is 

the  reduced  mass  between  species  a  and  r,  yg  is  the  ratio  of  specific 

heats  for  species  s,  i  is  the  mean  free  path  of  spe.  es  s,  and  Q  is  the 

8  >r 


collision  cross  section  between  species  s  and  r.  The  sums  are  over  all 

heavy  species.  Por  collisions  between  neutral  species,  the  collision  cross 

section  was  approximated  by  Q  _  -  np^  ,  where  p  is  the  Leonard -Jones 

sr  r  8  x  o 


parameter  between  speces  r  and  s.  Values  for  p  used  in  the  model  are 

Jf8 

listed  in  Table  XX-4.  Por  collisions  between  neutrals  and  ions,  a 

2 

collision  cross  section  of  10  A  was  used.  Por  collisions  between  ions, 
the  Spitzer  cross  section  for  charged  species  was  used. 


II  2 

Q  -  6fffZ  Z  1 
rs  l  r  sJ 
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where  A  is  as  given  in  Eq.  11-31. 


The  ratio  of  specific  heats  was  assumed  to  be  5/3  for  all  monoatomic 
species.  Por  the  diatomic  species  H2  and  H*,  7«c g/cw  was  obtained  by 
separately  calculating  cp  and  cy.  The  specific  heat  at  constant  pressure, 
c  was  obtained  from  the  formulas  and  constants  for  the  appropriate  species 
found  in  Ref.  XX. 5.  The  specific  heat  at  constant  volume,  c  is  calculated 

V 

by  noting  that  c  «  T(  3S/3T )  and  obtaining  the  entropy  S  from  the  formulas 


11—18 


and  constants  in  Ref.  il.s 
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TABLE  XX-1 


SPECIES  INCLUDED  IN  THE  MODEL 


Hydrogen: 


Argon: 


H? 

H2 

H(n),  n-1,10 


Ha<  v),  v-1.2 

H2#  €-8.85  eV,  12  ev 


Ar(n)  I,  n»l,ll 
Ar(n)  III,  n-1,6 
Ar(n)  V,  n-1,3 


Ar(n)  IX,  n— 1,7 
Ar(n)  IV,  n-1,5 
Ar  vi 


*V 


*V 
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Table  ii— 2 


ELECTRON  IMPACT  AMD  HEAVT  PARTICLE  COLLISION  RATES 


PROCESS 

RATE+ 

REFERENCE 

e  +  h2  -  h2  +  e 

a) 

1 

e  +  H2  -  H2(v-1,2)  +  e 

a) 

2,  b) 

*  ** 

e  +  H  -  H  '  +  e 

2  2 

a) 

3,  b) 

e  +  H2  -  2H  +  e 

a) 

4 

e  +  H  *•  H4  +  2e 

2  2 

€  «  15.4  ev  c),d) 

5,6 

e  +  h2*  -  2H  +  e 

e  -  12.0  ev 

•  ) 

e  +  H  *  -*  h!  +  2e 

2  2 

e  »  3.4  ev  c),d) 

5,6 

e  +  H4  -  H+  +  H  +  e 

a> 

7 

e  +  H(n)-*H(n)  +  e 

a) 

8 

e  +  H(n)  -  H(n+1)  +  e 

€n»13.6(l/n2-l/(n+l)2)  g) 

5,6,  b) 

e  +  H( n )  -  H+  +  2e 

2 

€  »  13.6/n  c),d ) 

5,6 

e  +  H4  -  2H 

e  +  Ar  -•  Ar  +  e 

a) 

10 

e  +  Ar(n)  -  Ar(n+1)  +  t 

9) 

5,6 

2  '  2+1 

e  +  Ar(n)  -  Ar(  1 )  +  2e 

5,6 

•  +  Xe  -  Xe  +  e 

a) 

10 

e  +  M2+  -  M2+  +  e 
f) 

2 . 9< -6  )2JP(  A  )/(  T2  /  2  )T  (  Z  ) 

11,12 

•  ♦  m2+  -  m<2-1>+ 
e  +  •  +  uZ+  -  m*  2-1  >+  +  e 
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H2+M-H  +  H+  M 

M  +  (H*,H2(v-1))  -  M  +  2H 
M+H2-M+H++H 
H2(v-0)  +  M  -  H2(v=1)  +  M 

M  +  H-M+  +  M+  e 

**  * 

H  +  M  -  H  4-  H 

**  * 

A r  +  M  -  Ar  +  M 


H=Ar  6.92(-6)T-1*°25esp(-51992/T  )  13 

M=H  2.95(-4)T93^2  3 

2  _3  /  2 

M=H  2 . 0(  —6  )T 

g 


to 


h) 


M*A r  4.153(-28 )T4*3 
M=H2  2 . 461 ( -27 )T* 

H»H  3 . 32(  -10 )exp( -1369/T  ) 


11 


l(-l6)vthexp(-€iATg)  i) 

2  >  j  \ 

H=Ar,H  1(-18)  can  J  ' 

M=H2  1( -17 )  cm 

2  \ 

M=Ar,H  1( -18  )  cm 
M=H2  1( —17  )  cm 


+  S(-7)  *  5  x  io~ 7  j  Rates  are  in  units  of  cm3-s-1  unless  noted  otherwise.  T  in  eV,  T 

e 

in  K 


a)  Cross  section  in  indicated  reference  was  averaged  over  a  Maxwellian 
electron  distribution  function  to  obtain  collision  rate  as  a  function  of 
electron  temperature. 

b)  Reverse  rates  are  by  detailed  balance. 

c)  Rate  is  a  function  of  threshold  energy,  €,  as  indicated. 

d)  Subject  to  lowering  of  ionization  potential. 

e)  Cross  section  and  rate  from  Refs.  5  and  6  were  used  where  the  threshold 
energy  is  equal  to  the  dissociation  energy. 

3  1/2 

f)  y(Z)  is  the  Spitzer  correction  term  and  A  -  1.55(10)(T  /n  )  '  .  For 
F(A),  see  Sec.  II-C. 

g)  For  threshold  energy,  see  Table  II-3. 

h)  Same  rate  as  for  H2  except  for  Boltzmann  factor  with  which  the  actual 
dissociation  energy  is  used. 

11-22 


/ 


References  for  Table  XX-1 


1.  D.  K.  Gibson,  Aust.  J.  Phys.  23,  683(1970) 

2.  G.  J.  Schulz,  Phys.  Rev.  135,  A983(1964) 

3.  A.  G.  Engelhardt  and  A.  V.  Phelps,  Phys.  Rev.  131,  2115(1963) 

4.  S.  J.  B.  Corrigan,  J.  Chem.  Phys.  43,  4381(1965) 

5.  H.  W.  Drawin,  European  CEA  Report  No.  PC-387,  1967;  European  CEA  Report 
No.  PC-383,  1966. 

6.  C.  Deutsch,  J.  Appl .  Phys.  44,  1142(1973) 

7.  C.  P.  Chan,  MReaction  Cross  sections  and  Rate  coefficients  Related  to 
the  Production  of  Positive  Hydrogen  Ions".  Lawrence  Berkeley  Laboratory 
LBID-632,  1983. 

8.  Y.  ItUcawa,  Planet.  Space  Sci.  19  993  (1971) 

9.  Erwin  and  Rung,  Trans.  Plasma  Science  PS-11,  266(1983) 

10.  L.  S.  Prost  and  A.  V.  Phelps,  Phys.  .  136,  A1538(1964) 

11.  M.  Mitchner  ar.-l  C.  H.  Kruger,  Jr.,  Partially  Ionized  Plasmas,  Wiley, 
New  York  (1973). 

12.  N.  Cohen  and  J.  P.  Bott,  "Review  of  Rate  Data  for  Reactions  of 
Interest  in  HP  and  DP  Lasers",  Aerospace  Corporation  fSD-TR-82-86,  Oct. 
1982 . 

13.  V.  N.  Kondratiev,  Rate  Constants  of  Gas  Phase  Reactions,  NTIS  COM-72- 
10014,  Jan  1972. 


TABLE  II-3 


Argon  Psuedo  state  Energies  and  Degeneracies 


Ar  I  Ar  II 


(I.P.=  126881. 

-1 

cm  ) 

(I.P.=  222434. 

-1 

cm  ) 

-1 

Energy  (cm  ) 

gtotal 

-1 

Energy  (cm  ) 

gtotal 

0. 

l 

0. 

4 

93450. 

8 

1432. 

2 

95000. 

4 

108722 

2 

104100. 

3 

132500. 

20 

105500. 

12 

135000. 

12 

106100. 

8 

139000. 

6 

107900. 

13 

157000. 

20 

111750. 

4 

113600. 

23 

116200. 

15 

119250. 

48 

Ar  nr 

(I.P.-  329303. 

-1 

cm  ) 

Ar  rv 

(I.P.~  481512. 

-1 

cm  ) 

-1 

Energy  (cm  ) 

gtotal 

-1 

Energy  (cm  ) 

gtotal 

0. 

5 

0. 

4 

14010. 

5 

21000. 

10 

33270. 

1 

34000. 

6 

115000. 

9 

117000. 

12 

144000 . 

3 

145000 . 

10 

187000. 

29 

\ 

s 
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Ar  V 

(I.P.-  604004. 


Table  II-3  (cont. ) 

cm  ) 


-l 

Energy  (cm  ) 


^total 


0.  1 

16300.  S 

121000.  15 

141000.  9 


Ref.:  C.E.  Moore  Atomic  Energy  Levels,  Vol,  I,  National  Bureau  of 
Standards,  NSRDS-NBS  35,  1971. 
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Table  II-4 

Lennard-Jones  Parameters  for  Calculation  of 
Collision  Cross  Sections 

Collision  Lennar d-Jone? 

Partners  Parameter  (A) 


A  -  A 


3.44 


A  -  H 


2.7 


A  - 


HL 


3.2 


2.97 


H 


2 


H 


2.6 


2.0 


Ref.:  J.o.  Hirschfelder,  C.P.  Curtiss,  and  R.Byron  Bird,  Molecular  Theory 
of  Gases  and  Liquids,  Wiley,  New  York,  1954. 
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SECTION  III 


SIMULATED  SPARK  COLUMN  CHARACTERISTICS 

In  this  section,  simulated  densities,  temperatures,  and  circuit 

quantities  for  a  spark  column  in  a  mixture  of  Xe/H2 1 0 . 01/0 . 9S  will  be 

examined.  The  initial  gas  p  assure  is  2  atm  and  the  charging  potential  is 

40  kv.  The  electrodes  are  separated  by  1.2  cm  and  the  current  return  path 

is  14  cm  in  diameter.  The  external  inductance  is  20  nH.  The  on  axis 

14  .3 

^ reionization  electron  density  is  10  cm  and  has  a  radius  of  50  Atm. 

The  duration  of  the  preionization  pulse  is  assumed  to  be  very  Short  so  that 
the  initial  electron  density  appears  instantaneously.  The  intent  of  this 
exercise  is  to  illustrate  the  canonical  characteristics  of  a  laser 
triggered  spark  column.  The  general  trends  discussed  here  apply  to  a  wide 
range  of  initial  conditions  and  gas  mixtures;  the  details,  of  course, 
change . 


III. A.  Hydrodynamic  Quantities 

The  density  of  heavy  particles,  in  unit's  of  &pfpo  are  is  plotted  in 
Fig.  Ill-la.  The  initial  density  is  2.7  x  10-  g/cm  .  The  same  data  is 
plotted  in  Pig.  Ill-lb  in  absolute  units.  The  pressure  and  and  gas 
temperature  fields  for  the  same  conditions  are  plotted  in  Figs.  III-2,  and 
III-3  respectively.  The  velocity  field  is  plotted  in  Fig.  III-4.  The 
radial  pressure  gradient  is  the  accelerating  force  for  the  convective 
velocity.  The  magnitude  of  the  accelerating  force  is  proportional  to  the 
density  of  contours  in  the  horizontal  direction  of  Fig.  III-Z.  The  gas 
remains  nearly  motionless  during  the  first  10  ns  as  this  is  the  time 
required  for  the  electron  density  to  avalanche  to  a  sufficiently  large 
value  to  begin  heating  the  heavy  particles.  Concurrent  with  the  increase 
in  gas  temperature  (Fig.  III-3)  the  gas  pressure  begins  to  increase,  first 
near  the  axis  and  later  at  points  at  larger  radii .  The  maximum  pressure 
within  the  arc  is  approximately  600  atm.  The  subsequent  pressure  gradient 
initiates  convective  motion  of  the  gas  iPig.  III-4). 

Convection,  and  subsequent  expansion  of  the  ionized  cere  of  the 
spark,  begins  as  the  gas  within  the  core  is  heated  by  electron  impact 
collisions.  Although  numerous  and  hot,  the  pressure  of  the  electron  fluid 
contributes  to  the  pressure  gradient  only  in  the  ratio  of  the  ion  mass  to 
the  total  heavy  particle  mass.  This  is  a  consequence  of  the  conditions  of 
ambipolar  diffusion.  (See  Section  II.  c, )  As  the  hot  gas  begins  to 
convect  radially  outward,  a  low  density  core  is  left  behind  and  a  region  of 
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increased  heavy  particle  density  begins  to  form  at  the  edge  of  the  ionized 

region.  (See  rig.  III-l)  Although  cold,  the  pressure  of  this  high  density 

gas  increases  to  the  point  that  it  too  begins  to  convect.  At  this  point,  a 

shock  begins  to  form.  The  shock  is  initially  strong  (A P/Pc  >5)  and  thin 

(width  <  5  tm),  and  later  weak  (A p/p  ■*  5 )  and  thick  (width  >  30  im). 

o 

The  mass  density  of  the  hot  core  is  reduced  to  a  value  <  0.0S  that  of  its 
initial  value,  and  the  core  is  nearly  isothermal.  The  inner  contour  in 
Pig.  XXI-3  corresponds  to  a  temperature  of  50,000  K.  The  temperature 
gradient  between  this  high  value  and  its  abmient  value  occurs  across  a  thin 
region  of  gas  just  interior  to  the  shock  and  penetrates  into  the  shock  for 
approximately  half  its  width.  These  conditions  impact  heavily  on  the 
distribution  of  atomic  and  molecular  species  in  the  plasma  column.  (See 
Sec.  XXX.  B. ) 

After  the  first  20  ns  of  convection,  the  pressure  gradient  within  the 
core  decreases  leaving  nearly  isobaric  conditions.  As  the  arc  continues  to 
expand,  pressure  in  the  isobaric  core  slowly  decreases.  Steep  pressure 
gradients  are  delegated  to  the  outer  regions  of  the  arc.  The  transition 
between  the  core  (pressures  of  many  hundreds  of  atmospheres)  and  the  cool 
gas  exterior  to  the  core  (pressure  of  a  few  atmospheres )  occurs  across  the 
shock  that  has  a  thickness  of  only  tens  of  microns.  The  location  of  the 
maximum  velocity  is  just  interior  to  the  shock  front.  At  a  particular 
radial  location,  the  convective  velocity  monotonically  increases  to  its 
maximum  value.  Although  remaining  positive,  as  the  convective  velocity 
decreases  from  its  maximum  value  it  may  oscillate,  responding  to  small 
changes  in  the  pressure  gradient.  These  pressure  and  velocity  oscillations 
result  in  oscillations  in  the  height  of  the  shock  front  (rig.  XXX-1),  the 
amplitude  of  which  damps  in  time. 

Convection  is  the  spark  column's  attempt  to  reduce  a  pressure 
gradient  The  convective  velocity  continues  to  increase  until  the  pressure 
gradient  is  reduced  to  a  low  value  as  a  result  of  the  redistribution  of 
mass  by  convection.  The  direction  of  the  convective  velocity  can  reverse 
with  a  sufficiently  large  change  in  the  sign  of  the  pressure  gradient. 

This  is  an  effect  which  is  commonly  seen  in  bounded  systems,  such  as  an  arc 
sustained  within  a  discharge  tube  where  convection  will  cause  mass  to 
accumulate  at  the  wall  thereby  increasing  the  gas  pressure  at  that 
location.  Xn  our  system,  though,  the  walls  are  sufficiently  far  away  that 
expansion  of  the  spark  column  during  the  first  100  ns  is  not  hindered  by 
significant  reversals  in  the  pressure  gradient.  Any  reversal  in  the 
pressure  gradient  is  not  large  enough  in  magnitude  or  long  enough  in 
duration  to  reverse  the  convective  velocity.  Under  certain  conditions. 


though,  significant  reversal  of  the  pressure  gradient  can  occur.  These 
conditions  will  be  discussed  below  in  Sec.  V.B. 

The  shell  of  high  neutral  particle  density  serves  two  nearly  opposing 
functions;  to  both  limit  the  rate  of  growth  of  the  column  and  to  provide 
the  seed  electrons  required  for  an  avalanche  the  gas,  thereby  expanding  the 
gets  by  non-convective  means.  The  accumulation  of  mass  in  the  shock  front 
increases  the  pressure  gradient  required  to  accelerate  that  mass  to  a  given 
velocity.  The  high  density  shell  also  presents  a  region  with  a  low  E/N 
which  inhibits  the  electrons  from  avalanching.  In  these  respects,  the  high 
density  shell  confines  the  arc.  The  high  density  shell  serves  as  a  source 
of  seed  electrons  by  trapping  ionizing  radiation  from  the  core  and  by 
thermal  ionization  resulting  from  its  adiabatic  compression.  The  fact  that 
the  photoionization  occurs  in  a  region  of  low  E/N  is  unfortunate  as  the 
probability  for  avalanche  is  low.  The  use  of  additives  having  an 
ionization  potential  below  that  of  hydrogen  and  xenon,  and  therefore  can  be 
ionized  by  radiation  to  which  the  shock  front  is  transparent  is  a  method  to 
increase  the  rate  of  arc  expansion. 


III. B.  Atomic  and  Molecular  Densities;  Ion  and  Neutral 

Electron  density  for  the  spark  column  examined  in  the  previous 
section  is  plotted  in  Fig.  III-5.  The  contours  span  the  range  from  2  x 
10  to  1  x  10  .  To  an  accuracy  of  a  few  percent  (see  below)  this 

Fig.  HI— 5  is  also  a  plot  of  the  of  the  density  of  H+.  The  electron 
density  reaches  its  maximum  value  near  the  axis  and  early  during  the 
discharge  pulse  before  the  heavy  particle  density  has  been  depleted  by 
convection.  The  maximum  electron  density  does  not  correspond  to  the 
maximum  in  fraction  ionized.  As  the  gas  convects,  thereby  lowering  the 
mass  density  in  the  core,  the  electron  density  decreases  by  a  smaller 
fraction  than  the  mass  density  thereby  increasing  the  ionized  fraction. 
After  35  ns,  the  core  is  fully  ionized.  For  gas  mixtures  that  are  greater 
than  95%  hydrogen,  the  electron  temperature  in  the  core  exceeds  15  eV. 

This  elevated  temperature  results  from  the  fact  that  after  the  molecular 
hydrogen  has  been  diriociated  and  the  atomic  hydrogen  ionized,  there  is  not 
a  significant  density  of  species  remaining  with  which  electrons  can  have 
inelastic  collisions.  The  only  remaining  energy  loss  mechanisms  for 
electrons  in  the  core  at  this  time  are  elastic  collisions  with  ions  or 
bremsthrallung,  neither  of  which  are  efficient  dissipators  of  energy. 
Cooling  of  the  electrons  occurs  as  the  core  expands,  lowering  the  impedance 
of  the  discharge  and  the  resulting  E/N. 


The  densities  of  xe+  and  X&++  (the  highest  ionization  state 
considered  for  this  calculation)  are  plotted  in  Pig.  III-6.  Xe  I  has  the 
lowest  ionization  potential  of  all  the  species  in  the  mixture  and  therefore 
is  the  first  to  be  ionized.  At  the  high  electron  temperature  afforded  by 
the  high  fraction  of  H2,  Xe  i  is  rapidly  ionized  to  create  a  region  where 
all  xenon  species  are  xe  XI,  Which  in  turn  are  ionized  to  leave  a  core 
where  virtually  all  xenon  is  Xe  ixx. 

Although  initially  the  most  abundant  neutral  species,  a2  is  also  a 
highly  volatile  species  in  the  hot  spark  column  environment.  The 
dissociation  potential  of  H2  is  4.5  eV  and  that  of  hJ  is  3.4  eV. 
Dissociation  proceeds  by  electron  impact  collisions  and  by  thermal 
dissociation,  requiring  only  moderate  heating  of  the  gas  (>  5000  K)  to 
result  in  significant  dissociation.  The  additional  mechansim  of 
dissociative  recombination  is  also  available  to  the  ion.  Once  the 
molecular  species  are  dissociated  a  major  electron  energy  loss  mechanism, 
that  of  vibrational  excitation  is  eliminated,  thereby  enabling  the  electron 
temperature  to  climb  to  a  higher  value.  The  higher  local  electron 
temperature  leads  to  rapid  and  ultimately  total  ionization.  This  sequence 
of  events  is  illustrated  in  Pig.  III-7  where  there  are  plots  of  the 
densities  of  S2>  H*,  and  H.  At  early  times,  before  the  gas  has 
significantly  heated  or  avalanched,  molecular  hydrogen  is  only  slowly 
depleted  near  the  axis.  After  about  20  ns,  molecular  hydrogen  begins  to 
first  ionize  and  then  dissociate.  To  the  left  of  the  low  density  contour, 
H2  is  fully  dissociated.  The  density  of  H2  is  high  in  the  shock  wave, 
however  the  width  of  the  region  in  which  H2  is  at  high  density  is  narrower 
than  that  of  the  shock.  (See  Pig.  XXX-lb.  )  This  indicates  that  H2  is 
thermally  dissociating  on  the  inner  side  of  the  shock  and  being  replaced  by 
atomic  hydrogen.  (See  Pig.  IIX-7b. ) 

After  the  first  30  ns,  r.tomic  hydrogen  exists  only  within  and  just 
interior  to  the  shock.  Xt  is  in  this  region  where  the  gas  is  hot  enough  to 
fully  dissociate  h2  however  the  electron  density  is  not  yet  high  enough  to 
deplete,  by  ionization,  the  atomic  hydrogen.  The  density  of  H*  proceeds  , 
through  a  similar,  but  even  more  volatile  sequence.  (See.  Pig.  III-?c. )1 
Early  in  the  discharge,  H*  is  the  most  abundant  ion.  As  the  gas  heats  and 
the  electron  density  increases,  H*  is  rapidly  dissociated.  The  region  ini 
which  the  rate  of  ionization  of  H2  exceeds  the  rate  of  dissociation  and  1 
recombination  of  H*  is  confined  to  a  very  narrow  portion  of  the  shock.  1 
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Figure  III-2  Thermodynamic  Pressure  for  a  Spark  Column 
in  a  Xe/H  :0. 01/0, 99  Gas  Mixture 


Figure  III-3  Heavy  Particle  Temperature  for  a  Spark  Column 
in  a  Xe/H: 0.01/0. 99  Gas  Mixture 


Figure  III-4  Velocity  Field  for  a  Spark  Column  in 
a  Xe/H  :0.01/0.99  Gas  Mixture 
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Figure  III-5  Electron  Density  for  a  Spark  Column  in 
Xe/H., :  0.01/0. 99  Gas  Mixture 


Figure  III-7  Density  of  Hydrogen  Species  for  a  Spark  Column  in 
Xe/H_:0.01/0.99  Gas  Mixture.  (a)  H  ;  (b)  H;  (c) 


SECTION  IV 


RATE  OP  ARC  FORMATION  AND  EXPANSION 

In  this  section,  the  rate  of  arc  formation  and  expansion  will  be 
examined  for  various  initial  conditions  including  preionization  electron 
density,  charging  voltage,  and  gas  composition. 

IV. A.  Delay  Time  for  Voltage  Collapse 

There  are  two  time  quantities  of  interest  in  laser  triggering  of 
spark  gaps.  The  first  is  delay  time,  the  time  between  the  laser  induced 
breakdown  and  voltage  collapse  of  (or  current  flow  through)  the  gap.  The 
second  quantitity  of  interest  is  jitter,  the  scatter  in  delay  time  for  a 
given  laser  intensity  for  a  large  number  of  trials.  The  former  quantitity, 
delay  time,  was  examined  with  the  model  as  a  function  of  initial  electron 
density.  These  results  are  plotted  in  Pig.  xv-1.  The  delay  time  for  this 
figure  is  defined  as  the  time  between  laser  triggering,  that  is  the  time  at 
which  the  preionization  laser  instanteously  produces  the  indicated  electron 
density,  and  when  the  current  reaches  a  value  of  100  A.  Clearly,  the  delay 
time  decreases  with  increasing  initial  electron  density.  True  avalancning, 
and  hence  current  flow,  does  not  occur  until  multistep  ionization  becomes  a 
significant  fraction  of  the  total  rate  of  ionization.  The  final 
characteristics  of  the  spark  column  »xe  relatively  independent  of  the 
initial  electron  density.  Obtaining  a  given  value  of  current  or  electron 
density  is  merely  shifted  in  time  by  a  value  equal  to  the  delay  time. 

To  compare  the  computed  values  for  delay  time  with  experimental 
measurements, an  assumption  must  be  made  as  to  the  laser  ionization 
process.  The  ionization  potential  of  the  neutral  gas  is  12. 5  eV  for  xe 
(mole  fraction  .01)  and  lb. 5  eV  for  H (mole  fraction  0.99).  The  energy  of 
the  laser  photon  is  **  4.4  eV  Therefore,  a  minimum  of  3  photons  are 
required  to  ionize  the  Xe  and  4  photons  are  required  to  ionize  the  B2<  If 
we  assume  the  ionization  to  be  a  '■-photon  process,  then  the  initial 
electron  density,  n^,  is  proportional  to  X  ,  where  X  is  the  intensity  of 
laser  used  for  preionization.  ^Dsing  this  assumption,  and  normalizing  an 
initial  electron  density  of  10  cm  to  an  intensity  of  10  W/cm  , 
experimental  points  are  also  plotted  in  Pig.  IV-1.  The  agreement  is  very 
good  indicating  that  the  ionization  process  may  indeed  be  by  4-photons. 
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IV. B.  Effect  of  Molecular  weight 


Recall  that  the  sound  speed  of  a  gas  at  temperature  T  and  of  average 
( rv  2  \ 

molecular  weight  M  is'  *  ' 


trv-i] 


Providing  that  the  dominant  expansion  mechanism  for  our  arc  columns  is  a 
hydrodynamic  sonic  expansion  of  the  hot  plasma,  the  expansion  rate  can  be 
effected  by  a  change  in  plasma  temperature  or  in  molecular  weight  of  the 
gas.  Arc  radius  (defined  as  the  outer  radius  at  which  the  electron  density 
falls  to  .05  of  the  instantaneous  maximum  value)  as  a  function  of  time  for 
three  different  Ar/H^  gas  mixtures  are  plotted  in  Pig.  IV-2  to  illustrate 
this  effect.  As  the  average  molecular  weight  of  the  mixture  increases,  and 
the  average  sound  speed  decreases,  the  rate  of  arc  expansion  decreases. 
Note,  though,  that  for  the  first  forty  nanoseconds,  which  includes  the  time 
delay  discussed  above,  the  rate  of  expansion  for  all  cases  is  nearly 
identical.  Althougn  the  heavy  particles  within  the  core  are  hot  (T^  > 

15000  K)  there  has  been  insufficient  time  for  the  pressure  gradient  to 
acclerate  the  particles  to  a  significant  velocity.  During  this  period, 
radiation  transport  and  subsequent  photoionization  of  the  neutral  gas 
outside  the  plasma  column  is  responsible  for  expansion  of  the  arc.  As  the 
gas  accelerates,  convection  becomes  the  dominant  expansion  mechanism.  At 
this  time,  the  expansion  rate  for  plasmas  with  different  me-  lecular  weights 
begins  to  diverge.  The  gas  mixture  with  the  lower  average  molecular 
weight,  and  corresponding  high  sound  speed,  expands  at  the  higher  rate. 


IV. C.  Recombination  as  a  Plasma  Heating  Mechanism 


Collirional  radiative  recombination  is  believed  to  be  the  dominant 
recombination  mechanism  for  the  arcs  of  interest.  The  rate  coefficinet  for 
this  process  is  inversely  proportional  to  the  9/2  power  of  the  electron 
temperature . ( 1  v * 3 )  Therefore  in  hot  plasmas  (many  electron  volts),  the 
recombination  rate  coefficient  is  small.  The  total  rate  of  recombination, 
though,  is  proportional  to  the  cube  of  the  electron  density.  Therefore 
even  in  hot  plasmas  with  large  electron  densities,  the  tctal  rate  of 
recombination  can  be  large .  in  a  collisional  radiative  recombination,  a 
portion  of  the  kinetic  and  potential  energy  of  the  recombining  electron  is 
carried  off  by  the  electron  that  is  its  collision  partner.  This  energy 
remains  in  the  electron  gas.  The  total  electron  density,  though,  is 
reduced  by  one  electron.  The  result  is  that  the  total  average  energy  of 
the  electron  distribution  increases .  This  process  is  known  as 
recombination  heating  and  is  a  mechanism  by  which  a  significant  amount  of 
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energy  is  retained  in  the  gas.  Recombinat ion  is  not  an  external  source  of 
energy,  as  is  joule  heating.  Rather,  recombination  heating  is  a  mechanism 
whereby  the  energy  of  the  recombining  electron  is  distributed  to  the 
remaining  electrons.  This  energy  is  originally  put  into  the  system  by 
external  means  and  is  stored  in  the  kinetic  and  potential  energy  of  the 
recombining  electron.  Had  the  recombination  been  purely  radiative  instead 
of  collisirnal  radiative,  then  the  energy  of  the  recombining  electron  could 
only  be  recovered  if  the  resulting  photon  is  absorbed  by  the  plasma. 

The  rate  of  recombination  heating,  R^,  is  given  by 

3 

*r  *  (lkTe  +  <1_f)€r  1  T/T 

T 

e 

Where  the  e_  is  the  ionization  potential  of  the  state  into  which  the 
electron  recombines  and  f  is  the  fraction  of  that  energy  that  is  radiated 
away  during  the  recombination.  The  first  term  of  Eq.  XV-2  is  for  the 
thermal  (kinetic)  energy  of  the  recombining  electron  and  the  second  therm 
is  for  the  potential  energy  of  th^  recombining  electron.  For  typical 
values  in  our  plasma  ( j  -  3  m/ cm  ,  a  »  100  O-cm)  the  rate  of  joule 
heating  is  90^/cra  .  The  rate  of  recombinat  ion  heating  (T^  -  4  eV,  » 

S  eV,  ne  «  10  )  is  approximately  IS  GW/cm  ,  and  is  an  appreciable 

fraction  of  the  total  energy  input.  Energy  by  this  process  can  have  a 
significant  impact  on  the  arc  dynamics. 

The  degree  to  which  recombination  heating  is  important  in  the 
expansion  of  the  arc  was  investigated  with  the  model,  the  results  of  which 
are  plotted  in  Pig.  IV-3a.  Arc  radius  is  plotted  as  a  function  of  time  for 
two  case  Which  differ  only  in  the  inclusion  or  absence  of  recomb inat ion 
heating  (Eq.  XV-2).  The  arcs  are  for  mixtures  of  Ar/H2t0.2/0.8.  Electron 
temperature  for  the  same  conditions  at  30  ns  after  triggering  is  plotted  in 
Fig.  XV- 3b.  The  rate  of  arc  expansion  xs  faster  for  the  case  with 
recomb inat ion  heating,  though  only  marginally  so.  The  electron  temperature 
for  this  case  is  only  about  10%  higher  near  the  axis  however  it  is 
siginifcantly  higher  towards  the  edge  of  the  spark  column.  (Note  the 
decrease  in  electron  temperature  in  high  density  region  where  E/N  is  low. ) 

Small  perturbations  in  the  temperature  of  the  plasma,  as  plotted  in 
Fig.  XV-3b,  can  cause  significant  changes  in  the  properties  of  the  plasma, 
especially  electron  density.  Recall  that  the  rate  of  ionization  scales  as 
exp(-cI/T#),  where  is  the  ionization  potential.  While  collisional 
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radiative  recombination  scales  as  T  .  Therefore  a  small  increase  m 

e 

the  electron  temperature  not  only  increases  the  rate  of  ionization  but  also 
decreases  the  rate  of  recombination.  Using  these  simple  scaling  laws,  one 
can  show  that 
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Electron  density  for  the  two  cases  discussed  above  is  plotted  in  Fig.  IV-4 
as  a  function  of  radius  approximately  30  ns  after  triggering  the  spark 
column.  The  case  with  recombination  heating  has  an  average  electron 
density  about  twice  that  of  the  case  without  reconciliation  heating,  and 
also  shows  more  "shell  structure"  than  the  case  without  heating.  This 
region  of  high  electron  density  corresponds  to  the  region  of  higher 
electron  temperature  in  Fig.  rv-3b.  (Note  that  both  cases  also  show  a 
penetration  of  the  plasma  into  the  high  density  shock  wave,  a  result  of 
local  sources  of  photo-  and  thermal-  idnization. ) 


§ 
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IV. D.  Arc  Expansion  as  a  Function  of  Charging  Voltage 

The  rate  of  arc  expans  ion  can  be  expected  to  be  a  function  of  both 
the  magnitude  of  energy  deposited  in  the  plasma  and  the  rate  of 
deposition.  Computed  tiass  density  and  arc  radii  for  spark  columns  having 
different  values  of  charging  voltage  are  plotted  in  Fig.  IV-5.  As  one 
would  expect,  the  rate  of  arc  expansion  increases  with  increasing  charging 
voltage.  Note  other  qualitative  trends  in  the  plots  of  mass  density  as 
charging  voltage,  and  hence  rate  of  energy  deposition  increaes.  The 
thickness  of  the  shock  decreases  and  its  height  (in  terms  of  A P/PQ) 
Increases  with  increasing  charging  voltage,  while  the  local  minimum  in 
shock  height  occurs  at.  progressively  earlier  times  with  this  increase.  The 
average  mass  density  with;'  n  the  core  of  the  spark  also  increases  as 
charging  voltage  increases.  This  trend  indicates  local  reversals  in  the 
pressure  gradient  resulting  from  rapid  heating  of  the  interior  of  the 
shock.  Reversal  of  the  pressure  gradient  can  slow,  or  reverse,  the 
radially  outward  directed  convection. 
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Figure  IV-1.  Computed  and  Experimental  Delay  Times  Retween 

Laser  Triggering  and  Voltage  Collapse  as  a  Function 


Figure  IV-2.  Arc  Radius  as  a  Function  of  Time  for  Different 
Molecular  Weights  and  Mixtures  of  Ar/H  . 


Figure  IV-3.  Conditions  With  and  Without  Recombination  Heating 
(a)  Arc  radius;  (b)  electron  temperature. 
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Figure  17-4.  Electron  Density  With  and  Without  Recombination  Heating 


Figure  IV-5  Mass  Density  (Ap/p  )  and  Arc  Radius  for  Spark 
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SECTION  V 


INVERSION  OF  COMPUTED  DATA  TO  SIMULATE  INTERFEROGRAMS 
AND  COMPARISON  WITH  EXPERIMENTAL  DATA. 

V. A.  Introduction 

Experimentally  obtained  inter ferograms  of  laser  triggered  spark 
columns  record  the  product  of  the  change  in  index  of  refraction  from  the 
ambient  value  with  the  path  length.  This  product  can  be  related  to  the 
local  electron  density  and  the  change  in  heavy  particle  density.  The 
measured  values  are  line  of  sight  (LOS)  data.  In  order  to  obtain  data  as  a 
function  of  radial  position  within  the  column,  an  Abell  inversion  mu3t  be 
performed.  For  LOS  data  recorded  as  a  function  of  the  transverse 
coordinate  x,  F(x),  the  corresponding  radial  function,  G(r)  is  obtained  by 

(V  n 

performing  the  integral  ‘  ' 


where  R  is  the  radius  of  the  cylinder  of  interest.  For  our  purposes,  the 
function  F(x)  is  t'..e  fringe  shift  of  the  probing  laser  resulting  from 
changes  in  heavy  particle  density  or  electrons.  In  order  to  separate  the 
contributions  of  electrons  and  heavy  part ic leu  in  the  inte ferograms,  and 
hence  obtain  their  densities,  inteferograms  at  two  different  wavelengths 
must  be  inspected,  and  their  inverted  profiles  subtracted.  In  principle, 
this  is  a  straight  forward  exercise.  In  practice,  though,  local  values  of 
the  fringe  shifts  change  on  time  scales  comensurate  with  the  duration  of 
the  laser  and  commensurate  with  the  jitter  (or  uncertainty)  in  measuring 
the  time  delay  between  initiation  of  the  spark  column  and  interferogram. 
Since  inverted  profiles  from  two  different  interferograms  taken  on 
different  spark  columns  (though  with  the  same  initial  conditioi-s)  are 
subtracted,  it  is  extremely  important  that  the  delay  time  of  the  two 
interferograms  be  the  same.  If  this  condition  is  not  met,  deconvolution  of 
the  fringe  shift  data  for  comparison  to  computed  results  is  difficult. 

An  equivalent  comparison  between  this  model  and  experimental  data  is 
to  construct  a  fringe  pattern  from  the  computed  densities  and  compare  the 
result  to  the  experimentally  obtained  interferogram.  This  process  requires 
one  to  compute  the  integral 
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where  now  F(x,t)  is  the  simulated  transverse  fringe  shift  at  time  t,  At  is 
the  time  duration  of  the  probing  laser,  X  is  the  wavelength  of  interest 
and  G(r,t)  is  the  change  in  index  of  refraction  from  unity  as  computed  with 
results  from  the  model.  G(r,t)  has  contribution  from  heavy  particles  and 
from  electrons.  For  heavy  particles. 
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where  p  is  the  mass  density,  PQ  is  its  initial  value,  and  A p  is  the 
change  in  mass  density  from  its  initial  value.  The  constant  y  depends 
weakly  on  the  gas  mix  and  for  our  purposes,  y  »  2  *  lo"  For 
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electrons,  is  given  by  the  plasma  dispersion  relationship*  *  '  ! 
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where  is  the  electron  density  and  in  cgs  units  &  »  8.95  x  lo  .i 

I 

The  integral  over  time  in  Eq.  V-2  is  necessary  because  during! the 
time  duration  of  the  laser  pulse  (5  ns),  the  shock  front  of  the  plasma 
column  can  move  a  distance  comparable  to  the  resolution  of  the  ! 
measurement.  Quantities  of  interest,  such  as  electron  density, $can j also 
change  appreciably  during  this  time.  With  a  velocity  of  5  x  10  cm/s,  a 
shock  wave  can  travel  a  distance  of  25  pm  during  the  laser  pulse.  The 
result  is  the  fringes  are  broadened  and  the  apparent  thickness  of  the  shock 
front  increased.  The  thickness  of  the  fringes  associated  with  the  j 
undisturbed  fringe  pattern  is  accounted  for  by  displacing  the  simulated 
fringe  patterns  by  a  distance  equal  to  the  undisturbed  thinkness. 


V.B.  Simulated  Xnterferograms  and  Comparison  to  Experimental  Data 


Experimental  inter ferograms  for  two  xe/H2  gas  mixes  are  in  Fig.  V-l. 
In  Fig.  V-la,  the  gas  mix  is  Xe/H2 : 0 . 01/0 . 99  and  in  Fig.  v-lb  the  gas  mix 
is  xe/H2: 0.03/0. 97.  The  charging  voltage  is  40  kV  and  gap  spacing  is  1.1 
cm.  The  undisturbed  fringe  pattern  is  at  the  extreme  right  or  left  of  the 


interferograms . 


Positive  displacements  axe  due  to  positive  A  p/p  ,  while 

o 
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negative  displacements  are  due  to  negative  Ap/p^  and  electrons.  Recall 
that  these  interferograms  are  of  line  of  sight  data.  The  fringe  shifts 
must  undergo  an  Abell  inversion  to  obtain  data  as  a  function  of  radius. 

The  different  appearance  of  the  fringes  for  a  relatively  small  change  in 
gas  mixture  is  dramatic.  In  the  LOS  data  of  Pig.  v-la  (Xe/H2: 0.01/0.99) 
there  appears  to  be  a  rather  wide  but  fairly  distinct  shell  of  high  heavy 
particle  density.  The  LOS  fringe  pattern  in  the  center  of  the  arc  appears 
somewhat  circular  and  is  characteristic  of  a  column  with  a  fairly  uniform 
electron  density  profile.  In  the  LOS  data  of  Fig.  V-lb  (Xe/H2s 0.03/0. 97), 
high  density  region  appears  less  distinct.  The  fringe  pattern  in  the 
center  of  the  arc  is  flatter  than  that  of  the  low  xenon  mixture  and  has 
approximately  twice  the  total  shift.  The  transition  to  the  center  plateau 
occurs  steeply  in  a  region  close  to  the  high  high  heavy  particle  shell. 

This  fringe  pattern  is  characteristic  of  “heavy"  gas  mixtures,  whereas 
"light"  gas  mixtures  tend  to  have  the  more  circular  fringe  pattern. 

Simulated  interferograms  using  results  from  the  SPRAD  code  for  the 
conditions  of  Fig.  V-l  are  plotted  in  Fig.  V-2.  The  qualitative  ag’T'ewnt 
is  excellent.  The  striking  difference  in  the  appearanc-  r>..'  the  fringe 
pattern  is  best  explained  by  referring  to  Fig.  v-3,  where  electron  density 
is  plotted  as  a  function  of  radius  for  the  two  cases.  For  the  low  xenon 
case,  the  electron  density,  confined  within  the  high  mass  density  shell,  is 
fairly  uniform  as  a  function  of  radius.  The  electron  density  penetrates 
into  the  shock  wave  with  only  a  moderate  slope.  In  contrast ,  the 
xe/H2 :0.03/0.97  case  has  a  high  electron  density  shell  interior  to  the 
shock  wave.  The  LOS  fringe  pattern  therefore  mimicks  that  of  a  thin  shell; 
a  flat  center  with  steep  sides.  The  shock  wave  of  the  Xe/H2: 0.01/0. 99  case 
is  somewhat  thinner  but  almost  a  factor  of  two  higher  than  that  of  the  the 
Xe/H2: 0.03/0. 97  case.  The  mass  density  in  the  interior  is  rather  uniform 
as  opposed  to  the  distribution  for  the  high  xenon  case  which  is  both  higher 
and  has  structure. 

The  differences  in  shape  of  the  electron  density  profile,  and  hence 
fringe  pattern,  between  the  1%  and  3%  xenon-hydrogen  mixtures  results  in 
put  from  the  difference  in  average  molecular  weight  (3.29  AMU  vs  5.87 
AMtJ),  and  therefore  sound  speed.  The  1%  xenon  mixture  is  sufficiently 
light  that  the  gas  responds  quickly  to  pressure  gradient.  As  a  result,  the 
convective  velocity  is  positive  ( radially  outward)  and  the  pressure 
gradient  is  negative  (pressure  decreasing  radially  outward).  The  molecular 
weight  of  the  3%  xenon  mixture,  though,  is  sufficiently  large,  and  the 
response  time  of  the  gas  sufficiently  long,  that  similar  conditions  cannot 
be  maintained .  ( See  Fig .  V-4 . )  As  a  result  of  the  lower  convective 

velocity,  the  high  density  region  interior  to  the  shock  is  heated  thereby 
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creating  a  local  maximum  in  pressure.  This  reverses  the  pressure  gradient, 
slows  the  convective  velocity  interior  to  the  local  maximum  and  may  infact 
reverse  the  particle  flux  to  he  radially  inward.  A  local  region  of  low 
particle  density  results.  This  region  has  a  correspondingly  higher  E/N  and 
therefore  higher  rate  of  ionization.  The  result  is  a  shell  of  high 
electron  density  interior  to  the  shock. 

For  certain  •'heavy'*  gas  mixtures,  especially  with  high  fractions  of 

xenon  and  late  in  time  during  the  discharge  pulse,  the  fringe  pattern 

exterior  to  the  shock  is  not  flat.  An  example  of  which  is  shown  in  Fig.  V- 

5.  Before  reaching  the  shock,  the  fringes  tilt  upward  in  the  direction 

indicating  Ap/pQ  >  1.  For  light  gas  mixtures  the  transition  at  the  outer 

edge  of  the  shock  is  abrupt;  that  is  the  fringes  remain  flat  until  the 

shock  is  reached.  A  simulated  interferogram  for  similar  discharge 

conditions  is  in  Fig.  V-6.  Results  from  the  model  indicate  that  the  tilted 

fringes  are  in  fact  a  result  of  having  A  p/p  >  1  exterior  to  the  shock. 

o 

For  the  lighter  gas  mixtures,  the  rate  of  energy  deposition  as  well  as  the 
temperature  excursion  across  the  shock  are  be-th  larger  thereby  resulting  in 
a  steeper,  thinner  shock.  For  the  heavier  gas  mixtures,  the  slower  rate  of 
energy  deposition,  as  well  as  a  smaller  total  fraction  of  deposited  energy 
being  available  for  the  thermodynamic  expansion,  results  in  an  apparently 
weakening  and  widening  of  the  shock  late  in  the  discharge  pulse. 
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Figure  V-3  Mass  Density  and  Electron  Density  for  Spark  Columns 

in  Heavy  (Xe/H,:0.03/0.97)  and  Eight  (Xe/H _:0.01/0.99) 
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Figure  V-6  Simulated  Interferogram  Showing  Anomalous  Fringe  Shift 
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SECTION  VI 


THERMODYNAMICS  OF  HIGH  TEMPERATURE  Ar/H^  GAS  MIXTURES  AND 
A  SIMPLE  MODEL  FOR  LASER  TRIGGERED  SPARK  GAPS 


VI. A.  Introduction 

From  the  results  of  our  detailed  modeling,  and  the  results  of  our 
companion  experimental  program^  ^ ' 1  \  we  have  learned  that  arc  growth  in  a 
laser  triggered  spark  gap  cam,  to  first  order,  be  described  am  a  convective 
expansion  of  the  hot  ionized  core  augmented  by  therms  1-  and  photo¬ 
ionization  at  its  periphery,  with  have  further  learned  that  the  interior 
of  the  arc  is  very  nearly  in  local  thermal  equilibrium,  with  these  facts, 
a  very  simple  model  of  a  laser  triggered  spaurk  gap  can  be  constructed. 

This  model  is  based  on  am  energy  balance  between  ioule  heating,  accumulated 
internal  energy  amd  energy  lost  by  radiation  during  the  duration  of  the 
spark.  To  construct  this  model,  the  thermodynamics  of  Ar/H2  gc.s  mixtures 
must  be  known.  Although  conceptually  only  a  solution  of  Saha's  equation, 
in  practice  obtaining  the  thermodynaunic  state  of  am  arbitrary  gas  mixture 
is  difficult.  A  computationally  simple  method  for  obtaining  this  solution 
is  discussed . 


Within  the  central  core  of  the  spaurk  column  of  a  laser  triggered 

spark  gap,  am  well  as  throughout  the  hot  region  of  conventional  arcs,  the 

thermodynamic  state  of  the  plasma  is  well  approximated  by  assuming  local 

thermodynamic  equilibrium  (LTE).  In  TTE  *'.e  electron  temperature  amd  the 

heavy  particle  temperature  atre  eqv*  —  densities  of  various  excited  or 

thermodynamic  states  of  a  paurticular  species  in  the  gas  axe  given  by  their 

/  yj  4  2  ) 

Saha  equilibrium  values.  For  charged  species, 
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where  n  is  the  electron  density,  m  is  the  electron  mass,  M  is  the 
density  of  species  M  in  charge  stati  Z,  and  F^(T)  is  the  paurtition  cunction 
for  species  M.  The  ionization  potential  of  the  higher  state  relative  to  the 
lower  state  is  Ej  and  the  lowering  value  for  the  ionization  potential  (see 
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below)  is  cL( n@, Z ) .  For  diatmcxnic  heavy  particle  species  that  may 
dissociate,  the  Saha  expression  for  partitioning  of  densities  is. 
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where  is  the  mass  of  the  monoatomic  species,  the  density  of  the  diatomic 
is  M2,  that  of  the  monoatomic  is  M,  and  the  partition  functions  are  FM2 
and  fm  respectively.  There  may  be  many  excited  states  permitted  within  a 
given  ionization  state  of  the  heavy  particle  species.  The  distribution  of 
these  excited  states  axe  given  by  the  Boltzmann  expression.  That  is, 

[Mij]  *  gijexp(l^ 


where  N. .  is  the  density  of  heavy  particles  of  species  i  in  excited  state 
j,  gi_  is  the  statistical  weight  of  and  £„  is  the  excitation 
potential  of  the  relative  to  the  ground  state  of  JT  .  In  principle,  the 
sums  in  Eq.  VI-5  diverge  because  there  is  an  infinite  number  of  excited 
states  as  the  dissociation  or  ionization  continuum  are  approached.  This 
problem  is  circumvented  using  two  methods.  For  species  in  Which  the  next 
higher  species  (in  energy)  corresponds  to  an  ion,  the  ionization  potential 
is  lowered  from  its  nominal  value.  (See  below.)  This  lowering  of  the 
ionization  potential  effectively  sets  an  upper  limit  to  the  quantum  number 
for  excited  states.  For  diatomic  species  in  which  the  summation  is  over 
vibrational  -  rotational  states,  the  summation  is  terminated  at  a  level 

commensurate  with  the  dissociation  energy  of  the  molecule.  (See  _ 

Sec.  VI. C. ) 

In  a  plasma  with  a  high  electron  density,  the  effective  ionization 
potential  of  atoms  is  lower  them  that  in  the  non-ionized  gem.  This 
lowering  of  the  ionization  potential  results  from  em  effective  broadening, 
and  hence  coalescing,  of  energy  levels  near  the  ionization  continuum.  Due 
to  the  difficulty  in  measuring  this  value,  only  theoretical  expressions  are 
available.  In  this  work,  the  lowering  value  for  the  ionization  potential, 
in  eV,  is  given  by 
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CJVZ)  -  6.96  x  10  n  '  S 


where  the  electron  density  is  expressed  in 


VI. B.  Method  of  Solution 


We  are  interested  in  obtaining  all  the  pertinent  heavy  particle 
densities  and  thermodynamic  properties  for  an  arbitrary  mix  of  Ar  and  H2 
when  specifying  the  total  mass  density  and  thermodynamic  temperature. 
Alternately,  we  would  like  to  obtain  the  same  quantities  when  specifying 
the  mass  density  and  internal  energy  of  the  plasma.  The  temperature  is 
then  a  variable  for  which  we  solve.  For  a  multi-component  plasma,  a  large 
system  of  non-linear  equations  consisting  of  the  appropriate  Saha  equations 
must  be  solved.  Traditionally,  large  matrix  inversions  are  performed  to 
achieve  this  end.  For  our  system  of  equations  where  we  specify  the  total 
mass  density  we  cam  exploit  the  principle  of  conservation  of  mass  and 
chaurge  in  order  to  obtain  an  iterative  solution  which  converges  rapidly. 

For  a  gas  with  a  total  density  of  Ar  atoms  of  Arm 
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where  Z  *  i  is  the  maximum  ionization  state  of  Ar  of  interest.  Mote 

DldX  I  a, 

that  the  Saha  equation  for  Ar  can  be  written  as 


.  ( 1-1 )+  ^  , 

Ari+  w  ^ _ D1  (T,ne) 
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where  D^(T,ne)  contains  the  ratio  of  paLrtition  functions  and  Boltzmann 
factor.  The  dependence  of  D|  on  electron  density  is  through  the  value  of 
ch  amd  therefore  is  weak.  Therefore,  in  an  interative  fashion,  we  can 
write 


n  d. 
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The  kemal  value  for  Ar  is  obtained  from 
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1  +  D, /n  +  D, D_/n  + 
1  e  1  2  e 


we  can  write  analogous  expressions  for  the  hydrogen  species.  For  a  total 
density  of  hydrogen  atoms  of  1^, 


H_  m  H  +  H+  2  •  H„  +  2  •  H* 
T  2  2 
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Note  that  Eq’s.  12  have  as  a  kernal  Hz.  Solving  for  the  value  of  H2  as  a 

function  of  h  ,  we  obtain 

T  2 

[Jr  || 

H2  “  2K  1  +  “2~  ‘I  *  P5[n.) 

■  J  J  j  . 
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1  /2 

where  J  **  F  '  (1  <■  Fj/ne>/2  and  K  »  (1  +  F^n^.  The  problem  is  then 
sol/ed  by  obtaining  the  kernals  from  Eq’s.  10  and  13,  which  reduces  to 
solving  for  the  electron  density.  This  last  value  is  obtained  by  noting 
that  for  our  system. 


n  -  H+  +  H*  +  £  iAr*+ 

6  2  i=l 
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The  iterative  solution  consists  of  specifying  ArT,  H^,,  and  an  initial  guess 

for  nfi.  New  values  for  the  ionic  species  are  ■  btained  by  first  obtaining 

values  for  the  kernals  from  Eq’s.  10  and  13,  and  a  new  value  for  n&  is 

obtained  from  Eq.  VI-14.  This  process  repeats  until  the  solution  converges 

to  the  desired  accuracy.  Only  tens  of  iterations  are  usually  required  for 

—  6 

this  method  to  converge  to  an  accuracy  of  10 


VI. C.  Partition  Functions 


The  partition  function  for  a  diatomic  molecular  species  is  given  by 
the  sum*''1-3* 
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Where  u,x,y,B,D,a,  and  0  are  the  usual  spectroscopic  constants 
eeeeee  e  + 

and  are  listed  in  Table  VI -l  for  H  and  H_.  The  values  of  J  and  v 

2  2  4.  max  max 

for  H2  were  obtained  from  Ref.  vi.4.  For  H2#  the  value  for  for  a 


given  v  was  obtained  from  solving 
B[v,Jmax]  '  «°'°>  '  Do 
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where  Dq  is  the  H2  dissociation  energy  (2.65  eV).  The  value  of  v^^  was 

obtained  by  setting  J  *  0  in  Eq.  VI-17  and  solving  for  v. 

max 

For  monoatomic  species,  the  partition  function  is 
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For  both  Ar  amd  H,  we  obtained  values  for  g^  amd  E ^  from  Ref.  VI. 5.  ^ For 
Ar,  we  calcuated  values  for  ionization  states  up  to  and  including  Ar  +. 

To  reduce  the  number  of  levels  to  a  managable  size,  excited  states  of  Ar 
were  grou; ,d  into  psuedo  states ^which  contained  all  levels  within  an  energy 
width  of  approximately  1500  cm”  .  The  statistical  weight  of  a  given 
pseudo  state  was  then  the  sum  of  the  statistical  weights  of  all  real  levels 
contained  within  the  energy  bin. 


VI. D.  Thermodynamic  Quantitities  amd  Transport  Coefficients 


After  solving  for  the  number  densities  of  the  different  species  in 
the  plasma,  various  thermodynamic  quantities  and  tramsport  coefficients  can 
be  obtained.  Three  quamtities  of  interest  aure  the  internal  energy, 
electrical  conductivity  and  radiation  energy  loss.  Assume  we  have  quasi¬ 
steady  state  conditions  and  that  out  time  scale  of  interest  is  long 
compared  to  the  thermal  equilibration  time  of  a  given  column  of  gas,  but 
short  compared  to  the  thermal  conduction  of  energy  exterior  to  the  column. 
For  these  conditions,  the  internal  energy  of  the  plasma  is  maintained  by  a 
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balance  between  power  input  to  the  column  turn  joule  heating  (proportional 
to  the  conductivity  of  the  plasma)  and  power  lost  from  the  column  by 
emission  of  radiation.  In  principle,  the  internal  energy  of  the  plasma  can 
be  arbitrarily  large  or  small  since  the  energy  balance  (in  the  steady 
state)  depends  only  on  the  t^uilibration  of  the  energy  input  and  radiation 
losses.  The  time  required  to  reach  the  steaay  state  value  of  the  internal 
energy  depends  on  the  initial  defect  in  the  rate  of  joule  heating  to  that 
of  radiation  loss. 


The  internal  energy  of  the  plasma  is  simply 
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where  is  the  density  of  species  i  excited  to  level  j,  is  its 
excitation  energy  relative  to  Ar( Z-0,  n*0)  or  H2(v-0),  and  y  is  the  ratio 
of  specific  heats.  The  terms  in  Eq.  VI-19  are  for  electronic,  vibrational, 
and  ionization  energy  stored  within  each  species,  translational  energy  of 
the  heavy  particles  and  translational  energy  of  the  electrons. 

To  obtain  the  rate  of  radiation  loss,  one  must  know  the  emissivity  of 
the  plasma.  For  a  cylindrical  plasma  of  diameter  d  and  length  L,  where  d/L 
«1,  the  emissivity  at  frequency  v,  c(v,d)  can  be  approximated  by^ VI  • 6 ) 

C(v.d)  -  1  -  expf— : - - - J  [VI-20] 

l  ll  +  o.7(a(  v)d)  '  +  o.08a(  v  )dJ  J 


where  the  photon  absorption  coefficient  is  <*(v).  With  €(v,d),  the 
volumetric  rate  of  energy  loss  from  the  cylinder  due  to  radiation,  Erad*  is 
obtained  by  integrating  the  Planck  function  over  wavelength; 
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Absorption  coefficients  used  for  this  calculation  are  identical  to  those 
described  in  section  II.  B. 


The  electrical  conductivity  of  a  plasma. 
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is  a  proportional  to  the  electron  density  and  inversely  proportional  to  the 
electron  collision  frequency.  The  electron  collision  frequency  is  the  sum 
of  two  terms,  collisions  with  neutrals,  v^,  and  collisions  with  ions, 
vcI .  As  the  ion  density,  and  hence  ion  collision  frequency,  becomes 
comparable  to  that  of  the  neutrals,  the  electrical  conductivity  asymptotes 
to  a  maximum  value,  called  Spitzer's  conductivity.  Once  this  limit  is 
reached,  from  ion  density  considerations,  the  electrical  conductivity  can 
only  increase  by  increasing  the  temperature  of  the  plasma.  For 
sufficiently  dense  plasmas,  called  non-Debye  plasmas,  the  assumptions  with 
which  the  ion  collision  frequency  is  derived  are  no  longer  valid.  The 
limit  of  validity  is  when  there  is  less  than  one  electron  within  a  Debye 
3phere,  that  is  the  mean  separation  between  electrons  is  less  than  a  Debye 
length.  When  this  limit  is  reached,  the  electrical  conductivity  decreases 
from  the  Spitzer  value.  In  the  Spitzer  formulation,  the  electron  collision 
frequency  in  a  plasma  with  temperature  T,  summed  over  ions  with  density  N. 
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The  last  factor  in  Eq.  VI-23,  InA,  is  obtained  assuming  that  A  *  1,  which 
is  the  condition  for  a  Debye  plasma.  The  first  order  correction  for  non- 
Debye  plasma,  that  is  when  the  inequality  involving  A  does  not  hold,  is  to 
replace  the  above  factor  with  its  exact  expression; 
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In  the  limit  that  A  »  1,  Eq.  VI-25  reduces  to  £nA. 

1 

J, 

VI. E.  Simple  Scaling  Model  for  Spark  Column 

Given  the  values  for  electrical  conductivity  and  rate  of  radiation 
loss,  a  simple  scaling  model  for  spark  columns  can  be  constructed.  For  a 
given  mass  density  and  gas  mix,  we  cam  calculate  oe(T)  and  E^T).  For  a 
total  current  I  conducted  through  the  spark  column,  having  current  density 
j,  and  diameter  d  (I  -  jnd  /4),  energy  conservation  requires  that 
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where  At  is  the  time  duration  of  the  spark.  Eq.  VI-26  states  that  the 
energy  dalferenre  re tween  joule  heating  and  that  which  is  radiated  away 
must  equal  the  internal  energy  of  the  plasma.  Recall  that  these 
relationships  critically  depend  on  the  cylindrical  nature  of  the  plasma  and 
the  validity  of  the  assumption  of  local  thermodynamic  equilibrium.  For  a 
given  total  current,  Erad(T,d),  Eint*T*  and  °e(T)'  E<5’  VI_2S  can  k® 
satisfied  by  a  set  of  values  for  arc  diameter,  d.  and  temperature,  T.  In 
order  to  obtain  a  unique  solution,  we  must  have  euiother  relationship 
between  d  and  T.  To  obtain  this  relationship,  a  model  must  be  constructed 
describing  the  expansion  of  the  arc  for  a  given  temperature.  Satisfying 
this  requirement  is  the  motivation  for  the  detailed  kinetics,  hydrodynamics 
and  radiation  transport  model  discussed  in  previous  sections.  For  purposes 
of  obtaining  such  a  relationship  to  apply  in  this  simple  model,  we  assume 
that  the  arc  diameter  results  from  a  combination  of  sonic  expansion  of  the 
initial  arc  column,  augmented  by  a  non-convective  component  due  to 
radiation  transport  and  thermal  shock  heating  of  neutral  gas.  Using  thin 
model,  the  final  arc  diameter  ar  a  function  of  temperature  of  the  plasma  is 
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where  0  is  the  non-convective  augmentation  of  the  arc  expansion  (equal  to 
1  for  a  purely  sonic  expansion),  d±  is  the  initial  arc  diameter,  y  is  the 
ratio  of  specific  heats,  and  M  is  the  mean  molecular  weight  of  the  gas 
mixture.  If  we  assume  that  as  a  result  of  the  non-convective  expansion  of 
the  arc,  mass  is  entrained  into  the  plasmr*.,  then  the  final  mass  density  of 
the  arc  (neccessary  for  the  calculation)  is 
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where  p  is  the  initial  mass  density. 

VI. F.  Results  from  the  Thermodynamic  Model  for  Mixtures  of  Ar/H2 

Number  densities  for  a  plasma  in  a  mixture  of  0.5/0.5:Ar/H2,  as 
calculated  with  the  thermodynamic  model,  are  plotted  in  Fig.  VI -1.  Tot 
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these  results,  the  mass  density  ^is  constant  ard  corresponds  to  1  atm  at 

room  temperature  (p  *  8.6  x  10  g/cm  ).  The  values  are  plotted  as  the 

ratio  of  number  density  divided  by  the  total  number  density  of  all  species 

in  the  plasma  including  electrons.  (Recall  that  the  total  number  density 

increases  with  increasing  temperature  as  a  result  of  ionization  and 

dissociation. )  As  the  temperature  increases,  the  density  of  Ar  I  decreases 

and  the  density  of  Ar  XI  increases.  As  the  temperature  increases  further, 

higher  ionization  states  become  the  dominant  Ar  species.  The  density  of  H2 

monotonically  decreases  as  it  is  thermally  dissociated.  The  density  of 

atomic  hydrogen  increases  at  low  temperatures,  a  result  of  dissociation  of 

H  ,  but  decreases  at  temperatures  greater  than  about  1  eV  due  to 
2  + 
ionization.  The  density  of  H2  is  never  very  large. 

The  electrical  conductivity  for  the  0.5/0.5tAr/H2  gas  mix  is  plotted 
as  a  function  of  temperature  in  Fig.  VI-2.  The  approach  to  Spitzer 
conductivity  is  clearly  indicated  as  the  neutral  molecular  and  atomic 
species  are  dissociated  and  ionized,  leaving  a  fully  ionized  plasma.  The 
Spitzer  limit  is  indicated  by  a  straight  line  on  the  log-linear  plot. 

Photon  absorption  coefficients  as  a  function  of  photon  energy  and 
temperature  are  plotted  in  Pig.  VI-3 .  The  atomic  levels  were  not  Stark  or 
Dopplar  broadened  in  calculating  these  values.  At  low  temperatures,  there 
is  structure  in  the  wavelength  dependence  of  the  absorption  coefficient  as 
a  result  of  ionization  of  excited  states  of  H  and  Ar  X.  At  the  higher 
tempertures,  we  have  two  local  maxima  in  aborption.  The  maximum  in 
absorption  at  low  photon  energies  is  due  to  inverse  Brshmst rallung .  For 
low  temperatures,  first  excited  ntaten  of  H  and  Ar  X  have  appreciable 
densities,  thereby  accounting  for  the  absorption  peak  around  4  eV.  As 
these  atoms  are  ionized,  and  inverse  Brehmstrallung  becomes  more  important, 
this  peal*  is  lost  to  the  continuum.  The  second  maximum  in  absorption, 
occuring  from  12  to  15  eV,  is  due  to  photo-ionization  of  ground  state  H  and 
Ar  X.  The  gradual  shift  of  the  threshold  value  of  this  peak  to  lower  values 
as  the  temperature  increases  is  due  to  the  lowering  of  the  ionization 
potential.  The  reduction  in  the  absolute  value  of  the  peak  at  higher 
temperatures  is  due  to  the  ionization,  and  subsequent  depletion  of  the 
ground  state. 

Power  radiated  from  the  plasma  cylinder,  for  USE  conditions,  as  a 
function  of  wavelength  and  temperture  is  plotted  in  Fig.  VI-4.  At 
sufficiently  high  temperature,  the  emission  appears  thermal  in  character. 
The  sudden  jump  in  power  loss  at  approximately  12  eV  is  a  result  of  the 
large  emissivity  for  photons  with  that  energy  and  above.  An  observer  far 
from  the  plasma  column,  measuring  the  photon  flux  through  a  depth  of 
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neutral  gas,  would  not  •  .*  these  photons  as  they  would  be  absorbed  in  the 
intervening  gas.  These  photons  are  absorbed  in  the  shock  front  which 
confines  the  expanding  plasma  column ,  and  are  in  fact  the  mechanism 
responsible  for  the  non-convective  component  of  the  expansion.  To  the 
extent  that  these  photons  are  initially  absorbed  in  a  region  of  gas  which 
does  not  carry  current,  the  UV  photons  can  be  considered  as  being  lost  from 
the  system  with  respect  to  the  energy  balance.  In  reality,  the  plasma 
col..  in  is  not  a  perfect  cylinder  but  has  a  gradient  in  temperature  and 
electron  density  between  the  plasma  core  and  cold  gas  exterior  to  the 
core.  These  UV  photons  are  absorbed  in  the  gradient  region  where  the 
current  density  is  small  as  compared  to  the  core.  Therefore,  the 
assumptions  made  for  this  analysis  remain  valid. 

As  discussed  above,  the  er?  rgy  balance  relationships  between  internal 
energy,  radiation  loss,  and  *  »e  heating  (Eq.  VI-28)  do  not  define  a 
unique  relationship  between  diameter  and  temperature  of  the  spark  column. 
The  expansion  model  (Eq.  VI-29)  is  required  to  obtain  this  relationship. 

An  example  of  these  relationships  are  in  Pig.  VI-5,  where  arc  radius  as  a 
function  of  diameter  of  the  plasma  column  is  plotted  for  total  currents  of 
1-10  kA.  The  gas  mix  is  0.5/0.5tAr/H2.  The  initial  pressure  of  the  cold 
gas  is  2  atm  and  the  initial  diasieter  of  the  spark  is  100  Jim.  The  non- 
convective  (ie.<  radiative)  augmentation  to  the  convective  expansion  is  0 
-  '.s,  and  the  spark  duration  At  »  125  ns.  For  a  given  temperature, 
larger  diameter  sparks  are  required  to  carry  larger  currents.  Analogously, 
at  a  constant  diameter,  a  higher  temperature  plasma  will  carry  higher 
currents  because  the  conductivity  increases  with  increasing  temperature. 

T.  ere  is  a  large  motiviation,  with  respect  to  minimizing  the  loss  of  energy 
from  the  plasma,  to  maximize  the  diameter  of  the  column.  The  rate  of 
energy  loss  due  to  radiation  is  proportional  to  the  surface  area  of  the 
column  while  the  rate  of  joule  heating  is  proportinal  to  the  volume  of  the 
column.  Therefore,  having  a  plasma  column  with  a  larger  diameter  reduces 
tM  relative  loss  of  energy  to  radiation.  This  enables  a  higher 
temperature,  large?  conductivity,  and  higher  current  carrying  ability  for  a 
given  amount  of  joule  heating. 

Pig.  VI-5  also  contains  a  plot  of  Eq.  VI-29,  the  model  for  arc 
expansion,  which  defines  the  relationship  between  diameter  and 
temperature.  The  intersection  of  the  line  of  this  plot  with  those  lines 
defining  the  allowed  values  of  d  and  T  (based  on  energy  conservation)  for  a 
given  value  of  current  establishes  the  plasma  conditions  for  that 
particular  current.  Arc  diameters  and  temperatures  for  different  mixtures 
of  Ar/H2,  obtained  in  this  manner,  are  plotted  in  Pig.  VI-6.  The  initial 
pressure  is  2.0  atm.  As  the  mole  fraction  of  B2  increases,  the  arc 
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diameter  increases  as  a  result  of  the  reduced  average  molecular  weight 
having  a  correspondingly  larger  sound  speed.  In  the  Spitzer  limit,  the 
conductivity  depends  only  on  plasma  temperature.  Therefore  to  carry  a 
given  current ,  an  arc  of  smaller  diameter  ( ie . s  large  mole  fractions  of 
argon)  must  be  at  a  higher  temperature.  Analogously,  as  the  mole  fraction 
of  hydrogen  increases  and  diameter  of  the  arc  increases  (due  to  the  larger 
sound  speed),  for  a  given  total  current  the  required  plasma  temperature 
decreases . 

Arc  diameters  and  plasma  temperatures  with  an  initial  mix  of 
0 . 5/0 . 5 ! Ar/H^  as  a  function  of  current  are  plotted  ir.  Fig.  VI-7.  Results 
for  initial  pressures  of  1-3  atm  are  plotted.  Since  the  speed  of  sound 
depends  only  on  the  ratio  of  specific  heats  and  molecular  weight,  and  not 
on  density,  the  diameter  of  the  arc  changes  little  as  a  function  of 
pressure  for  a  given  current.  What  change  there  is  due  primarily  to  a 
change  in  temperature.  For  a  given  current  the  plasma  temperature 
increases  with  .^creasing  pressure.  When  the  intial  number  density 
decrease-,  the  ion  density  corresponding  to  the  Spitzer  limit  represents  a 
larger  fractional  ionization  thereby  requiring  a  temperature.  To  obtain 
this  fraction,  a  higher  temperature  is  required. 
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TABLE  VI-1 


Spectroscopic  Constants  of  H2  and  H2 


_i 

(cm  ) 


w 


w_x 
e  e 


cj  y 
exe 


8(J  x 
e  e 

a  u 
e  e 

u 

B 

"  3 

e 

e 

24Be 

♦395. 2 

117.99 

0.29 

60.809 

2.993 

36,045. 


/?  -  D 

e  e 


D  ■»  D„  +  ^  CJ  -  \  <J  x  +  i  U  v 
e  o  2  e  4  e  e  8  e  • 


2297. 

62. 


29.8 


1.4 


21,324. 


t  Ref.j  G.  H.  Herzberg,  Molecular  Spectra  and  Molecular  Structure 
I.  Spectra  of  Diatomic  Molecules,  2nd  Edition,  Van  Nostrand,  New  YorX, 
1950. 
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Figure  Vl-2.  Electrical  Conductivity  for  an  Ar/H2 : 0 . 5/0 . 5  Gas 
Mixture  as  a  Function  of  Temperature. 


PHOTON  ENERGY  (eV) 


Figure  VI-3  Optical  Absorption  Coefficient  for  an 

Ar/H2:0. 5/0.5  Gas  Mixture  as  a  Function 
of  Temperature  and  Photon  Energy. 
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Figure  VI-4.  Radiative  Energy  Loss  from  a  Spark  Column  in  an 
Ar/f^tO.S/O.S  Gas  Mixture  as  a  Function  of 
Temperature  and  Photon  Energy. 
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Model  Results  for  a  Spark  Column  for  an  Initial  Pressure  of 
n  Different  Ar/H2  Gas  Mixtures  as  a  Function  of  Total  Current 
c  diameter;  (b)  arc  temperature. 
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